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Commissioning of the new shock tunnel T4 at the University of Queensland 
implied that it was no longer necessary to focus the work of the research 
group about an annual test series conducted in the T3 shock tunnel in 
Canberra. Reflecting this, it has been possible to organize the group for 
work to proceed along lines such that particular personnel are associated with 
particular project areas. 

Thus the format of this report consists of a series of reports on specific 
project areas, with a brief general introduction commenting on each report. 
The introduction is structured by project areas, with the title of the 
relevant report stated under the project area heading. The reports themselves 
follow in the order of the project area headings. 
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1. COMMISSIONING OF SHOCK TUNNEL T4 


"The University of Queensland Free Piston Shock Tunnel T4 - Initial Operation 
and Preliminary Calibration" 

During 1987, the free piston shock tunnel T4 was brought into operation, and a 
preliminary calibration was performed. This involved selection of a 

particular driver gas volumetric compression ratio (X » 60) , and measurement 
of shock speeds and pressure history at the downstream end of the shock tube 
after shock reflection. These measurements yielded the stagnation enthalpy 
and pressure in the reservoir of test gas, at the end of the shock tube, which 
was supplied to the hypersonic nozzle. 

Preliminary test section calibration was effected through pitot rake 
measurements and static pressure measurement on a flat plate. Since the 

nozzle was a replica of one which has been used for some time with the T3 
facility, the results of pitot surveys taken in that facility were used. The 
ratio between the static pressure and the nozzle stagnation pressure then 
could be used to confirm the accuracy of the chemical non-equilibrium nozzle 
flow code used to predict the test section conditions, with an effective 
nozzle area ratio obtained from the ratio of pitot pressure to nozzle 
stagnation pressure. 

Since the static pressure is sensitive to helium contamination, the variation 
of static pressure with time could also be used as a check on contamination 
free test times. It was found that the measured times exceeded those 
predicted. 

With this calibration, it was cosidered that the tunnel was ready for some 
preliminary studies of hypersonic combustion. 


2. HYPERSONIC COMBUSTION STUDIES 

"Hypersonic Combustion of Hydrogen in a Shock Tunnel" 

Although static pressure levels were limited to approximately 20 kPa, it was 
found that it was possible to produce hypersonic combustion in a duct of 25 mm 
x 51 mm cross section, with a length of 600 mm. 
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As might be expected in a hypersonic flow, boundary layers appear to play a 
significant role. Not only do they generate high local temperatures, and 
thereby possibly influence ignition temperatures, but they also lead to 
significant pressure rises within the duct. Evidence of the first was 
apparent in the low ignition temperatures experienced in the experiments, and 
evidence of the second was apparent in measured pressure rise along the duct 
in the absence of combustion. Since the Reynolds' numbers of the experiments 
are such as to raise the possibility of the existence of turbulent, or transi- 
tional, boundary layers, it is important to learn more about transition under 
the conditions used in the tunnel. 

It is possible that the production of hydrogen combustion in the duct at Mach 
numbers from 4.8 to 6.3 was due to these boundary layer effects. 

When the duct was operated with the downstream part of one wall divergent, in 
order to produce thrust, it was found that specific impulse values up to 800 
sec were produced. This value occurred at a stagnation enthalpy of 
approximately 11 MJ/kg, confirming that hypersonic combustion would allow 
higher specific impulse at high stagnation enthalpies than supersonic 
combustion would allow. It was found that, unlike supersonic combustion, the 
peak specific impulse occurred at a stagnation enthalpy well in excess of the 
ignition value. Noting that the flow conditions achieved by the tunnel 
appeared to be somewhat marginal for hypersonic combustion, it is difficult to 
advance reasons for the slow increase in thrust with increase in temperature, 
since it could be influenced by reaction kinetics, boundary layer pressure 
gradients and other effects. However, it seems clear that hypersonic 
combustion offers new phenomena which must be understood. 


3. SIDEWALL ITJTL INJECTION 
Experiments 

"Hydrogen Scram jet with Sidewall Injection. Shock Tunnel Simulations" 

In previous work, it has been shown that injection of hydrogen fuel at the 
sidewall of a combustion duct is very effective in cooling the wall, but that 
observable pressure rises due to combustion did not occur until fuel 
equivalence ratios approached values of 2. 
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The present study is an examination of results of experiments performed in an 
produce thrust with a wall injected fuel layer. For this purpose, 
two stage nozzle expansions were tried, anticipating that a preliminary, low 
angle expansion would allow combustion to develop, whilst simultaneously 
developing thrust, and this thrust would be added to that produced by the 
second stage of expansion. However, only small .improvements were observed over 
the rather low levels of thrust produced by a single stage of expansion. 

The experimental results suggested that the cool walls of the model may be 
lowering the temperature of the mixing gases, and thereby would inhibit 

combustion. A numerical study was undertaken to explore this matter further. 

Numerical Study 

"Numerical Modelling of Sidewall Injected Scramjet" 

A numerical model, based on the CHABNAL computer code, was used to calculate 
the flow field for two conditions under which tests were conducted with a 
constant area duct, and one in which a single stage nozzle expansion was 

included. The code employed the k-c turbulence model, and an 8 reaction model 
for the non-equilibrium process of hydrogen combustion. The validity of the 
model was confirmed by comparing predicted values of surface heat transfer and 
pressure distributions with experimental results. 

The numerical simulation showed that wall quenching was not a significant 

factor in limiting combustion. In fact, the reason for the weak combustion 
effects observed in the experiments is the limited mixing which occurs between 
the hydrogen and the air streams. For the two cases studied (i.e. stagnation 
enthalpies of 4.2 MJ/kg and 8.7 MJ/kg respectively), the degree of mixing 
achieved 400 mm downstream of injection represented only 10% and 25% 

respectively of the injected hydrogen flow. The hydrogen which did mix with 
the air was almost completely burned. 

In simulations with injection of heated hydrogen, the amount of mixing 
increased in the lower enthalpy case, but not for the higher enthalpy. The 
combustion induced pressure rise was increased in the first case, but reduced 
in the second case, corresponding to changes in degree of mixing. 

Temperature contours in the hot hydrogen jet show that the hot regions of the 
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flow extend almost the entire way to the wall, indicating virtually no effect 
due to wall quenching. Thus the numerical simulation indicates that the 
temperature of the wall is not a significant factor in the experiments. 

It might be noted that in the experimental studies, the substantially 
increased pressure rises observed when a 20:80 Silane-Hydrogen mixture was 
injected at a stagnation enthalpy of 4.2 MJ/kg were interpreted as evidence 
that it was the depression of temperatures, rather than mixing limitations, 
which was limiting the combustion of hydrogen alone. Given the low level of 
mixing evident in the numerical studies, it now seems plausible that the lower 
injection velocities of the higher molecular weight Silane-Hydrogen mixture 
could have led to an increase in mixing which was sufficient to produce the 
observed pressure rise. 

The simulation suggests that mixing might be improved, and combustion effects 
increased, by using turbulence generators on the surface upstream of the 
injector. In fact, the moderate success achieved, in fig. 9 of the experimental 
paper, with transverse jets in that location could be interpreted as being due 
to such an effect. 


4. CHEMICAL KINETICS - EFFECT OF fREE STREAM ATOMIC OXYGEN 

"Combustion of Hydrogen and Hydrogen/Silane Mixtures" 

This study involves the use of one dimensional flow models, in which chemical 
reaction models are incorporated in a relatively complete form, in an attempt 
to explain effects observed in previous experiments in terms of chemical 
kinetics, and in assessing the effects of free stream atomic oxygen. 

For combustion of pure hydrogen, the effect of free stream oxygen concen- 
tration on ignition delay lengths was most pronounced at low temperatures, 
with the presence of atomic oxygen concentration, measured as only a 
fractional percentage, theoretically causing considerable reaction to take 
place in a flow where it otherwise would not be observed. As temperatures 
increased above 1000 K, the differences in ignition delay lengths reduced, and 
at high temperatures, even the presence of 2% atomic oxygen did not 
significantly alter the distance required to complete combustion. Clearly, the 
presence of free stream atomic oxygen has negligible influence on combustion 
lengths only at relatively high intake temperatures. 
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For si lane /hydrogen mixtures, previous chemical kinetics mechanisms were 
re-examined with respect to the role of silane in promoting ignition. Silane 
does this by increasing the temperature, and by increasing free radical 
concentrations. At temperatures below 800 K, scavenging of H0 2 also is 
important . The role of free radical concentrations suggests that atomic 
oxygen may have a less important role to play than in the case of hydrogen. 
Indeed, this is borne out by the computations, with the influence of free 
stream atomic oxygen being reduced, but nevertheless, the same trends are 
evident as with hydrogen, with the greatest sensitivity to atomic oxygen 
occurring at the lowest static temperatures. 


The one dimensional models produced results which were approximatley in accord 
with experiments, except at intake static temperatures of approximately 500 K. 
There ignition distances obtained experimentally are much shorter than those 
predicted, whilst the pressure rise is not as large. Also, higher pressures 
were achieved with lower silane concentrations. These results could not be 
explained, although it was pointed out that the limited pressure rises 
experienced suggested a mechanism involving only partial burning of the fuel. 


5. SCALING STOP US 
Experiments 

"Pressure Scaling Effects in a Scramjet Combustion Chamber" 

A re-examination was made of the results of experiments with a simple scramjet 
model in which the effects on net specific impulse of the pressure in the 
model, and the length of the combustion duct, were studied. 

With a given combustion duct length, the pressure dependence effects could be 
classified into three regions. In the first region, temperatures were too low 
for combustion, whilst in the second region, vigorous combustion took place, 
and the level of specific impulse was strongly dependent on the pressure. The 
boundary between the two regions occurred at about the same temperature, for 
pressures in excess of one atmosphere. As pressures were reduced below that 
level, the ignition boundary moved to higher temperatures and, as pressures 
approached 25 kPa, no combustion was observed at any temperature. In the 

third, higher temperature, region, specific impulse was not strongly dependent 
on pressure level. 
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In varying the duct length, it was found that a short duct at a high pressure 
could produce the same specific impulse as a longer duct at a lower, 
appropriate, pressure. However, when the pressure in the longer duct was 
raised to the high pressure level of the shorter one, it was observed that 
specific impulse was the same for the two ducts over most of the range of 
temperatures studies. It was concluded that, at the higher pressure level, the 
flow was approaching complete combustion in the duct. 

A numerical model give results which were reasonably consistent with the 
experiments in region 3, but not in region 2. However, the model was useful 
for indicating trends in this region. 

The results were consistent with a model of combustion which allowed for a 
pressure sensitive region, in which the ignition and reaction lengths were a 
dominant chemical feature of the flow, and a region which approached chemical 
equilibrium. However, because Mach number and Reynolds' number also varied as 
pressures were varied, it was not possible to be more specific about the 
factors influencing the flow. This indicated a need for an experiment to be 
designed in which the variables likely to affect combustion could be varied 
independently. Also, because of the ultimate importance of geometrical sealing 
if laboratory experiments are to be applied to flight situations, the 
experiment should also include this factor as a variable. 


Dealgn of a Scaling Experiment 

"Pressure Scaling in the Scramjet Model" 

Since the highest thrust levels had been achieved at temperatures just above 
ignition values, it seemed likely that scaling effects in this range would be 
of most interest. Therefore, in order to bias the experiments towards this 
objective, it is planned to ensure that the product of pressure and a typical 
model dimension is held constant. Three models are to be made, of sizes such 
that the geometrical scaling factor varies by a factor of 10. The experiments 
will initially involve variation of temperature for each of the three models, 
aiming- to obtain data for which Mach number, free stream composition and 
Reynolds' number are held constant between the three geometrically similar but 
different scale models. 
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In order to do this, it was necessary to use the same nozzle reservoir 
conditions and the same nozzle throat size for all three models. This ensured 
that free stream composition remained constant but, to ensure that the Mach 
number remained the same for the three models for differing static pressure 
levels, it was necessary for two of the models to expand the test flow to Mach 
numbers m excess of the combustion value, before compressing it through 
oblique shocks as it passed into the combustion duct. 


The experiments can be repeated at different Reynolds' numbers by changing the 
shock tunnel operating pressure levels. In order to change the Mach number, it 
is necessary to make a new shock tunnel nozzle to operate with the smallest 
model, but the other two models allow adjustment of Mach number through 
adjustment of the angle of the oblique shock generators. 


These experiments, and associated analysis, are expected to proceed over the 
next two or three years. 


*. PC? AMS IOW TUBE S TOPICS 

Experiments 

"Experiments on an Expansion Tube with a Free Piston Driver" 

Experiments were performed on a small expansion tube which employed a free 
piston driver, in order to determine if the flexibility offered by a free 
piston driver would allow the range of operation of an expansion tube to be 
extended. An existing free piston driver was used, and an expansion tube was 
constructed which was to be a representative V« scale model of the Langley 
expansion tube, with a free piston driver attached. Thus, both facilities 
exhibit a combined shock tube and acceleration tube length which is 
approximately 140 times the internal diameter. However, since discussions with 
Langley personnel indicated that there may be some advantage in using a longer 
shock tube than for the Langley facility, the ratio of shock tube length to 
acceleration tube length was 0.65 for the University of Queensland facility, 
compared with a value of 0.27 for the Langley one. 
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In the experiments, it has been shown that 

(i) The pitot pressure records on the expansion tube centreline are 
qualitatively similar for the Langley facility and the Queensland 
facility, when both are operated at similar Reynolds' number. However, 
at the higher acceleration tube filling pressures, an irregularity 
occurs in the Langley results which is not evident in the University of 
Queensland results. 

If this irregularity is ignored, then the results indicate that the 
University of Queensland facility is providing a reasonable simulation 
of the behaviour of the Langley facility. 

(ii) Test times of at least 80 fcisec were obtained, during which the 
fluctuations in pitot pressure were limited to ±5% or less. It was 
found that the test section velocity at which such test times could be 
obtained could be changed by changing the driver gas. Thus, a test 
section velocity of approximately 9.5 km/ sec was obtained with helium 
driver gas, and 4.5 km/ sec with argon driver gas. 

(iii) The overall operating pressure levels were reduced by a factor of 2 
{i.e. main diaphragm burst pressure, shock tube filling pressure, and 
acceleration tube filling pressure were all reduced by the same factor) 
whilst driver gas compression ratio was held constant and it was found 
that the effects on the pitot pressure tracer of varying acceleration 
tube pressure were qualitatively similar to those at the higher 
pressure levels. In particular, a useful test time was produced at the 
same pressure ratios. Further tests were done with the operating 
pressure levels reduced by another factor of 2, with the same results. 

These results indicate that overall operating pressure levels can be 
decreased or increased by a substantial factor without changing test 
times . 


9 



10 . 


Theory 

"Expansion Tube Test Time Predictions” 

Theoretical studies were made of the influence of the driver gas-test gas 
interface in the shock tube on the test time in the acceleration tube. Two 
mechanisms were considered for producing the pitot pressure fluctuations 
which limit the test time. The first involves the formation of low density 
gas inclusions ("blobs") within the test gas, due to driver gas-test gas 
mixing, and acceleration of these inclusions into the test region through the 
unsteady test gas expansion in the acceleration tube. The second involves 
reflection of the unsteady expansion from the turbulent interface region 
between the driver gas and the test gas . 

An unsteady method of characteristics code was written for a McIntosh desk 

computer, and was used to predict the arrival of each of the two types of 
disturbances in the test region. For the NASA Langley expansion tube, and the 
University of Queensland expansion tube with argon driver gas, it was 
predicted that the reflected expansion disturbances would be evident on the 
test section pitot trace before the inclusion disturbances. In both these 
cases experiments showed disturbances arriving somewhat in advance of the 
predicted time but nevertheless, showing a trend which was compatible with 

predictions. For the University of Queensland expansion tube with helium 
driver gas, the same effects were evident for the lower shock tube filling 
pressures used. However, for the higher shock tube filling pressures, the 
theory predicted that the inclusion disturbances would arrive first at the 
test sections. This was supported by the experiments, which showed 

disturbances of larger amplitude than the reflected expansion ones, arriving 
at the approximate time predicted. 

The encouraging results obtained with this preliminary theoretical analysis 

suggest that it has an important role to play in explaining the performance 
limitations of expansion tubes, and that it is worthwhile extending the 

analysis in order to test it further against experiments. 
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1 . INTRODUCTION 


The University of Queensland free piston shock tunnel T-4 is the latest in a 
series of free piston shock tunnels, the first three of which were constructed 
at the Australian National University, Canberra, in the 1960's. The complete 
series is shown to scale in fig. 1, and is marked by an increase in sire in 
passing from one to the next. 

The increase in size in passing from T3 to T4 was associated with proposed 
studies in hypersonic combustion and propulsion. In order to produce static 
pressures in a hypersonic flow which are high enough to allow vigorous 
combustion of hydrogen to take place in the relatively short lengths available 
in a shock tunnel model, it is necessary to produce high nozzle stagnation 
pressures. This means that these high pressures must be produced at the 
downstream end of the shock tube, and sustained for long enough to establish 
and maintain a test flow in the nozzle for an adequate period. 

A disappointing feature of T3 was the loss in pressure which occurred in the 
shock tube.* 1 * Although the cause of this loss is not understood, it has been 
found that it can be overcome by using a long compression tube for the free 
piston compressor. Hence the piston travels approximately 26 m during the 
driver compression stroke in T4, compared with 6 m for the T3 piston. 

T4 began operation in April, 1987 and, after a settling in period, came into 
routine operation in September, 1987 . At the time of writing, some 200 test 
runs of the facility have been made, and main diaphragm operating burst 
pressures have been raised to 28% of the design value. Experiments have shown 
that this is sufficient to allow combustion of hydrogen at a Mach number of 5. 
Although further increases in operating pressure levels are planned (in order 
to produce more vigorous hydrogen combustion) , it seems appropriate to use the 
test data obtained so far for’- a preliminary review of the operating 
characteristics of the facility. 
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2 . EXPERIMENTAL ARRANGEMENT 


As noted above, the free piston driver is 26 metres long, and it is 230 mm in 
diameter. The piston mass is 92 kg. The shock tube is 10 m long, and 76 mm in 
diameter, with shock timing stations 3, 2 and 1 located 2.00 m, 4.00 m and 
6.00 m respectively from the downstream end of the tube. The shock speed 
between station 2 and 3 was used for performance calculations, as it was 
considered that the gas processed in that part of the tube would be the test 
gas in the test section during the steady flow period. 

A contoured nozzle was used, with a throat diameter of 25 mm and an area ratio 
of 100. This produced a test section Mach number which varied from 4.8 at a 
stagnation enthalpy of 35 MJ/kg to 6.3 at 3.7 MJ/kg. The change in Mach 
number is due to the real gas effects in the nozzle expansion. 

3. SHOCK TUBE PERFORMANCE 


For shock tunnel purposes, the important aspects of the shock tube performance 
are the stagnation enthalpy and the stagnation pressure of the test gas 
supplied to the nozzle, as well as the time between shock reflection and 
arrival of driver gas contamination at the entrance to the nozzle. The 
latter, of course, determines the effective shock tunnel test time. 

The results reported here apply to one test condition, in which the 
compression tube was operated to produce a helium driver gas volumetric 
compression ratio of 60 (based upon the compression tube filling pressure to 
the main diaphragm burst pressure, and the assumption of isentropic 
compression) , and the main diaphragm burst pressure was 57 MPa. 

The nozzle stagnation pressure was obtained from two piezoelectric pressure 
transducers (PCb piezotronics type 118A) which were located in the shock tube 
wall, 65 mm upstream of the downstream end of the shock tube. The stagnation 
enthalpy was calculated by using the measured shock speed and the shock tube 
filling pressure to obtain the enthalpy and pressure immediately after shock 
reflection, and then assuming an isentropic expansion to the measured 
stagnation pressure. 

Test results are presented in fig. 2, with a scale for the measured shock speed 
matched to the shock tube diaphragm pressure ratio at the left hand side of the 
figure. Noting that Earth orbital velocity corresponds to a stagnation 
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enthalpy of 31 MJ/kg, it can be seen that stagnation enthalpies in excess of 
this value have been achieved in the shock tube. This is in accord with 
^< 2 > 

experience in T3 

The nozzle stagnation pressure was measured 0.5 milliseconds after shock 
reflection and, after being normalized with respect to the main diaphragm 
rupture pressure, has been plotted at the top of fig. 2. The tendency for this 
pressure to increase with falling stagnation enthalpy is consistent with 
expectations based on coupling shock tube theory with the observation that 
tailored interface operation occurs at a shock speed of 6.65 km/ sec., corres- 
ponding to a stagnation enthalpy of 40 MJ/kg. 

The variation of the nozzle stagnation pressure with time is shown in fig. 
3(a) for typical test conditions. For fig. 3(a) it can be seen that the 
pressure falls by approximately . . % over a period of one millisecond after 

shock reflection. It is thought that this is due to the use of a piston with 
a mass which was designed for operation at high pressures. At the pressures 
of tests, the piston velocity was too low at diaphragm rupture to match the 
flow of driver gas into the shock tube, causing the driver pressure to fall. 
This can be rectified by lowering the speed of sound in the driver gas, as in 
fig. 3(a) (ii) . Here it can be seen that, when argon is used as driver gas, 
rather than helium, the pressure remains approximately steady for a much 
longer period. 

Fig. 3(b) shows the limitations on test time arising from contamination of the 

test gas by the driver gas. This can come from two causes. One is due to 

drainage of the test gas through the nozzle throat, and the other is due to 

" jetting" of the driver gas along the walls of the shock tube from the contact 

surface. The "jetting” is due to the bifurcated shock pattern which forms when 

( 3 ) 

the strong reflected shock interacts with the shock tube boundary layer 
It can be see that the shock boundary layer interaction is responsible for 
limiting the tunnel test time. It will be noted that the test time has been 
expressed as an approximate length of the slug of test gas passing through the 
test section by multiplying the calculated test time by V2H^ . The test 
slug lengths which are displayed are based upon a conservative calculation and, 
as will be seen below, there is evidence to suggest that test times may be 
somewhat longer than those calculated. 
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4. TEST SECTION CAL I BRAT I OK 

Pitot survey data are shown in fig. 4. These were obtained in shock tunnel T3, 
using a nozzle which was identical with the one employed in these tests. They 
show that the nozzle expansion is not complete for some 0 inches downstream ‘of 
the expected nozzle station, yielding pitot pressures which are somewhat above 
the uniform flow values in a region which is downstream of the theoretical test 
cone boundary. This may be due to boundary layer effects in the nozzle flow, 
which were not taken into account in designing the nozzle. Downstream of this 
region, the pitot traverses indicate a uniform and parallel flow. 

Static pressures were measured at a pressure orifice located 200 mm from the 
leading edge of a flat plate mounted in the test section with its leading edge 
at the nozzle exit. A static pressure record, normalized with respect to the 
nbzzle stagnation pressure, is shown at the top of fig. 5. 

Static pressure measurement provides a good indication of the static 
temperature at the test section, since the velocity is essentially fixed by the 
stagnation enthalpy, and the density therefore is fixed by mass flow 
considerations. The temperature in the test section is sensitive to real gas 
effects in the nozzle expansion, and therefore it follows that the static 
pressure measurement serves as a check on the validity of the numerical 
calculations used to predict the test section conditions. In fig. 5, static 
pressure measurements are seen to compare satisfactorily with prediction made 
using a one-dimensional non-equilibrium calculation^ 

For comparison, the static pressure predicted for a perfect gas expansion of 
air to the same pitot pressure level also is shown in fig. 5, and is seen to be 
well removed from the measured values. Expansion of helium to the same pitot 
pressure levels would produce an even lower pressure. In this light, the decay 
in static pressure after about 1 millisecond of steady pressure, which is 
evident in the static pressure record in fig. 5, may be seen as the onset of 
driver gas contamination. This would indicate that the steady test flow slug 
lengths indicated in fig. 3 may be conservative by as much as a factor of two. 
However, this should be regarded as a tentative conclusion, which is subject to 
confirmation when more data is‘ obtained. 

It may be noted that some points are included on fig. 5 which were taken in the 
intake duct of a scramjet model. The static pressure there was higher than the 
measured flat plate static pressure but, assuming that the proportionality 
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between flat plate and intake static pressures would apply at all stagnation 
enthalpies, the intake results were scaled to produce psuedo flat plate 
measurements at stagnation enthalpies where flat plate measurements were not 
available 

The consistency between measured and calculated static pressures indicates that 
the numerical code may be used to calculate other test section flow properties. 
Results of such calculations are presented in fig. 6. The Mach number decays 
as stagnation enthalpy is increased because of real gas effects in the nozzle 
expansion. Two curves are shown for the free stream enthalpy. This is because 
a substantial part of the free stream enthalpy, h_, is invested in the 
sensible enthalpy associated with the temperature of the gas and, if desired, 
this could be converted to freestream kinetic energy by further expansion. 
However, the rest of the freestream enthalpy which is denoted by h f , is 
associated with "freezing" of the air as it passes through the nozzle 
expansion, and cannot be recovered by further expansion. It can be seen that 
this component rises to . 20% of the stagnation enthalpy at high stagnation 
enthalpies . 

5. CONCLUSION 

These results and calculations are preliminary in nature, representing an 
initial calibration of the tunnel at one range of operating conditions. 
Increases in operating pressure levels and changes in driver gas volumetric 
compression ratio are planned to generate new test conditions, as well as more 
extensive calibration measurements. However, the results reported here 
indicate that the shock tunnel is operating satisfactorily, at pressures and 
stagnation enthalpies which have been found to be high enough to allow 
experimentation on hypersonic combustion of hydrogen. 
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Hypersonic Combustion of Hydrogen in ft S tl PCk T unncl-L. 

Introduction. - 

This paper reports the results of a preliminary series of tests 
in the shock tunnel, T4, on hypersonic combustion and thrust 
generation in a hydrogen scram jet model. The work was performed 
before a full calibration of the tunnel and nozzle was completed, 
but using a nozzle profile . which had previously been tested in 
the shock tunnel T3. The main purpose of the experiments was to 
confirm the feasibility of hypersonic combustion m the flow 
conditions created by the new tunnel. 

Previous tests with a supersonic nozzle performed in the Shock 
tunnel T3 at the Australian National University, Ref 1, indicate 
that severe performance degradation occurs at the enthalpies 
associated with high Mach number flight. This appears to be at 
least partly due to the high static temperatures created by 
maintaining a supersonic combustion chamber at hypersonic flight 
speeds. Good combustion efficiency will require low combustion 
chamber intake temperatures, less than 2000 K, and this can only 
be achieved at enthalpies higher than 13 Mj/kg, corresponding to 
a flight speed of approximately 5km/s, by maintaining hypersonic 
flow in the combustion chamber. Even at lower enthalpies it is 
likely that hypersonic combustion will be preferable to 
supersonic if it can be made to burn properly. 

The methodology of the experiments was to initially perform a 
short series of tests in a constant area duct, to see if 
combustion was possible at all with the apparatus. The duct had 
an intake section with a height of 25 mm and a width of 51 mm . 
Injection was from the centreline across the full width of the 
duct. A thrust surface was then fitted at a suitable distance 
down stream of injection and specific impulse measurements were 
obtained for a range of conditions in the mach number range 4.8 
to 6.3. 

Compression wedges were then fitted to the model intake, and the 
flow in the combustion chamber was reduced to supersonic. With 
conditions adjusted to maintain the same intake pressure, a 
direct comparison was then obtained between supersonic and 
hypersonic combustion. The supersonic combustion tests were 
within the nominal mach number range of 4 to 4.5. 

After the completion of these runs a short series of tests was 
performed using an enlarged constant area duct with a nominally 
square intake section of 51 by 49.5 mm. The purpose of this was 
to distance the mixing layer from the presence of the wall 
boundary layer in an attempt to enhance combustion. 


Constant area d uct, hypersonic combustion. 

A schematic of the model used for these tests is shown in Fig 1. 
It incorporates the same injector strut as had been previously 
used with a supersonic nozzle, and its leading edge was extended 
clear of the intake so that no disturbances entered the duct. 
Injection was through a nozzle with a 1.6 mm throat on the 
centreline extending across the full duct width. Excepting 
disturbances from the side wall, the flow was "considered to be 
two dimensional. Instrumentation consisted of a series of 
pressure transducers on the top and bottom walls, both upstream 
and downstream of injection. 


25 



freMtp!l* m V,aS + in '^ ec !' ed from a room temperature reservoir with 

£3? 1 £ c Sr pr ks 

.£7 . ® ,, 5 nas a continuous pressure droD with ti». 
this is reflected in the nozzle .v<t *?„/, 2 p • h time, and 

rlnlilti shook? Ui as "given sS 

?ha 5 ^r? d t° ?* ' St '" dy ’ "hen^e ra?io P of sHtic p«sJi?e ll 

he nozzle to stagnation pressure, measured in the reservoir ic 
operating* 1 at t !! iS m j ght indloat « that the nozzle is 

lo ^at^oz^rLl^ rlTe^ii^ 

tSr.t^r e oi e<i the Th ^r?^r ? ressu L e 

ps-s 

pressure «eord a?e"°sho^ . Sf*S*T‘° n 

STl,^ rir Wh ? n H th ' *:L^Sre°^'aS 

flow is presented in Sa^oiM^f^cSterbM^TO^th^tiB^whSn 
the normalised pressure traces are steady. me when 

SiH£” iyi V e p“issi^ 

the°same ^n^W^" t0 ^luTS 

^d S lf"7 °Mj/S bUS as°Sa? r be e seerfrom “?J 1Pi “g°2 p ^!5 dV 

substantially S?SS ?or the' "oth^^’ tat - thd 

ST'and 1 ? le'd r<1 ? Se J hS dUCt i hteke r mach S 'unb« varies 'between 

rates. As th«e° "o^tioL^orJ^pc^o^^nlSe"^? t°" 
temperatures of between 2500 and 3600 K i-m* static 

reduced combustion efficient due m Lw ' due t L the 
dissociation in this temperature rjge “islZr. 

ssa^rf^ , «as 

Additionally, the intake velo^” i^ Mgh ii°!SI i22! r 5 t 2^!- 

mixing^^the”?;:? ^1““? 5«s dU 1t t [r S ^ S r jw.^at^m^ri^ 


in this enthalpy range, 
between approximately 6.1 , 

combustion chamber mach numbers, 
not will depend on the static 
proposed flight altitude. 


which corresponds to‘ flight speeds of 
and 7.3 km/s, will require even higher 
Whether this may be achieved or 
pressure levels relevant to the 


Further tests were performed in the enthalpy range of 13 to 6 ? 

££• iSSTSS-t^* 1 *^ x;? 1 **' v b ~ 3 d ^ 

5.81. In Fig 3 it can be seen that pr^SSME^ ZtSucZiZX 

26 



3 

injecting hydrogen fuel, over and above that produced by the 
injection of helium. In Fig 4 the fuel on traces are compared for 
a range of different enthalpies. The enthalpy range was not 
extended to identify the low temperature ignition limit, but it 
is noted that at the lowest enthalpy used for the constant area 
duct, 6.2 Mj/kg, the intake static temperature was down to 765 K. 
At this temperature, and at the intake pressures used, hydrogen 
fuel would not be expected to react within the confines of the 
duct without some sort of ignition enhancement. 

This is possibly provided in' a hypersonic duct by boundary layer 
heating. In Fig 8 the temperature profiles produced by 
equilibrium flow at mach 5 for laminar and turbulent boundary 
layers are shown. The flow conditions chosen correspond to 
the 8.39 Mj//kg stagnation enthalpy condition. In both cases the 
boundary layer was allowed to develop for a distance of 160 
mm, which corresponds to the upstream length of the injector 
strut. Whilst equilibrium calculations cannot be expected to 
exactly represent the state of a boundary layer in non 
equilibrium flow, it does illustrate that temperatures of the 
order of twice the free stream values are to be expected. The 
presence of this hot layer of gas might make it possible to 
sustain combustion at free stream conditions where it would not 
otherwise be possible. Furthermore, the presence of these hot 
spots might lead to local dissociation of the oxygen molecules, 
providing a source of radicals to encourage ignition. 

Another possible explanation is suggested by Fig 6, where the 
fuel off normalised static pressure distribution along the duct 
is shown, along with the predicted intake level. At an enthalpy 
of 8.39 Mj/kg the pressure is seen to agree well with the 
predicted value for a short distance after the injection strut, 
but further downstream the measured pressures are of the order of 
twice the theoretical values. This might be due to boundary layer 
growth in a restricted duct compressing the freestream flow. 
Pressure levels are uneven due to waves reflecting across the 
duct and from the walls, but the increased overall pressures 
downstream may clearly be seen. At this mach number, with the 
presence of waves indicated by the wall static pressure profile, 
it is hard to tell accurately what the local static pressure 
level will be in the free stream and in the mixing layer away 
from the wall. However the wall pressures are increased over such 
an extended region that it must be concluded that increased 
pressures also occur in the mixing zone. At mach 5, assuming any 
pressure increase over intake conditions occurs isentropically, 
a rise of approximately 20% will be produced in the static 
temperature when the static pressure doubles. 

Also shown in Fig 6 is the fuel off pressure distribution for an 
enthalpy of 26.5 Mj/kg. It can be seen that the downstream 
pressure does not consistently rise in this instance, but 
fluctuates around the predicted value. Assuming the pressure rise 
seen at 8.39 Mj/kg is indeed induced by boundary layer growth, 
then this would be qualitatively consistent with the reduced 
displacement thickness associated with laminar boundary layers at 
the lower reynolds number, higher enthalpy condition. It is 
possible that transition is occurring at the higher Reynolds 
numbers. Simple calculations may be used to indicate if the 
pressure rises observed are consistent with laminar or turbulent 
boundary layer growth by considering the compressive effect of 
the displacement thickness in the duct. 
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of the duct and the wave processes indicated by the pressure 
profiles are unknown. The calculations were repeated for a 
laminar boundary layer at the higher enthalpy of 26.5 Mj/kg. 
Although the lower Reynolds number produces thicker boundary 
layers, the displacement thickness does not increase 
significantly due to the higher ratio of static temperatures 
across the boundary layer. The physical reason for this may be 
seen from Fig 11 where the non dimensional density profiles 
across the boundary layer are shown. As both conditions share the 
same wall temperature of 3O0K, the cooled gas near the wall will 
have a higher non dimensionalised density in the high enthalpy 
flow, due to the high free stream static temperatures. This will 
serve to reduce the mass flow deficit of the boundary layer, and 
therefore the displacement thickness does not increase as would 
be expected in incompressible flow. 

Laminar boundary layers at the higher enthalpies, with transition 
to turbulence occurring at the higher Reynolds numbers associated 
with the lower enthalpy conditions is a possibility. However, 
even the presence of a laminar boundary layer at the higher 
enthalpy conditions would be expected to produce a pressure rise 
of the order of 303S in the duct, and as can be seen from Fig 2 
this was not observed. This might be because the approximate 
boundary layer model used assumes perfect gas behavior, with 
representative properties chosen at an appropriate mean 
temperature. This would lead to overpredicting the temperature, 
and underpredicting the density, in the boundary layer 'hot 
spots' where viscous heating predominates over heat loss to the 
wall. This in turn would lead to overestimating the displacement 
thickness by an amount which would be expected to increase with 
enthalpy. 

As a summary of this section, hypersonic combustion was achieved 
in a constant area duct for flows with stagnation enthalpies 
below 13 Mj/kg. At enthalpies above 13 Mj/kg no measurable 
combustion effects were observed. The high temperature cut off 
point is somewhere between 13 and 18.7 Mj/kg. The low temperature 
cut off point has' not yet been identified. An ignition 
enhancement mechanism is apparent, which is provisionally being 
attributed to the presence of boundary layers on the injector 
strut and on the side walls. Three modes are suggested whereby 
boundary layers might stimulate ignition. Firstly, the boundary 
layers on the injector contain regions of elevated temperature 
in close proximity to the fuel jet which might directly 
accelerate ignition and reaction. Secondly, the boundary layers 
may contain dissociated oxygen radicals which have to reach a 
certain concentration before the reaction 'takes off' . These two 
mechanisms are inherently interrelated. Finally, the growth of 
all boundary layers in the duct presents a restriction to the 
flow of the main stream and causes both pressure and temperature 
to rise. A preliminary look at the measured pressure profiles 
indicates that turbulent boundary layers may have been present 
for some of the conditions where combustion was observed. 


Comparison of Hypersonic and supersonic combustion With — & 

diverging duct. 

A configuration with a short injector with a leading edge located 
downstream of the intake was chosen for the main series of tests 
of hypersonic combustion in the thrust producing mode. The short 
injector was used in the hope that the extra intake area made 
available by having the leading edge compression process internal 
would increase the duct pressure enough to improve combustion. 

29 


6 


geometries p ^ ofiles a « plotted for the two 

pressure*^ i-h^ l bh * lnternal injector does show higher 
pressure in the constant area section, it is not 

downstream of the corner. The internal strut produces Uniform 

flow SUr ThA fte t the expansion, as would be expected in supersonic 

corner -Pol To strut shows a smaller pressure drop round the 

corner, followed by a gradually falling pressure downstream TKic- 

une ven 1 f 1 ovT^f i ^ numbe r ' with belong injector and a more 

f ii; p ° ssibly to the extra boundary layer 

growth m the middle of the duct. Despite an apparent reduction 
in performance the totally enclosed configuration was used for 

sat i sf Tot lndCr ° f th ? i ests as i-fc appeared to give a more 
satisfactory supersonic flow. 

In Fig 15 the two fuel on profiles are shown for equivalence 

ratios of 0.55 and 0.57 for internal and external leading edges 

in^ector Ve i y mH lt . oan be seen that th ® long strut of the extended 
tntTTw - Produced more pressure rise and thrust than the short 
totally enclosed injector. This may again be due to the larger 

edtfe d of y th? /er enbrained lnto the mixing flayer from the trailing 
i t the injector. Despite the extra mass flow induced with 
the internal injector, the pressure levels in the early part of 

at thi? 1 MSch a ^ er K d ° n Sj- re ^ lect this due to the wav e trajectory 
wh - + - cb nu ®ber. This is illustrated schematically in Fig 16 

not'ref lect “iVtTtE* 'T* ^ th * injector edfe di 

after injection. h Centre llne untl1 approximately 130 mm 

fSr FiS th« 3 “l- 14 t T C fUe 3 on . *** off Pressure profiles are shown 

SSS iM-S s ^ 

In Figs 9 and 10 the results are presented for hypersonic * n H 
created" by • , T !! e sup ' rsonic test conditions were 

7 5 dirges The di.ch^' S °°? P SiVine a double reflection of 

• a " grees - lhe diaphragm rupture pressure was *d inef .h 

aa ?i caPd V°t* ma 7he y ibdector^^as^configured 0 "^^^ "^th^^am^ 
hypersonic'JtsSlts^TablS l ff shlws%h dl nlSle C ^i? r cSSditions fS? 

bbe . .^^P® rs °o ic tests. Table 3 shows the supersonic test 
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for hypersonic and supersonic ducts respectively. 

Looking firstly at Fig 10, specific impulse against static intake 
temperature, it is seen that there is little difference between 
hyper and super sonic combustion. The low temperature cut off is 
less clearly defined in hypersonic flow, possibly due to more 
pronounced boundary layer heating. This effect would be hard to 
separate from the reduced particle residence times at the higher 
mach number. From a qualitative point of view it might be said 
that as the temperature cools down the higher mach number flow 

u S heat “ ip bas less time to and react within 

the duct However, as it does contain a small amount of very hot 

gas in the boundary layer, a residual amount of fuel still gets 
burned even at very low temperatures, but it does not release 
enough heat to ignite the main body of the fuel. 
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The very high levels of specific impulse previously observed at 
Mach 3.5 just above the low temperature ignition limit were not 
seen here. Noting that previous tests at Mach 3.5 showed the 
size of this peak to be very pressure sensitive, the lack of high 
specific impulse is probably due to the low static pressures in 
the hypersonic duct. It should be noted that combustion in this 
regime of potentially high specific impulse, that is between 400 
and 1500 K, cannot be strictly described as either a diffusion 
flame, where reaction times are considered insignificant, or as 
a premixed flame where heat release is totally reaction 
controlled. Overall heat release is therefore sensitive to any of 
the parameters effecting ignition, reaction and mixing, and is 
not fully understood at present. However, speculation as to 
the physical mechanisms dominant may be useful. 

The presence of a wall alongside the combustion zone of a 
scramjet with a combustion mach number of 3.5 has been shown to 
seriously reduce the effectiveness of combustion. The reason for 
this is not known, but it is possible that in a hypersonic flow, 
with thicker boundary layers, the same wall interaction effect 
may propagate into the mixing layer of a central jet too. If this 
is the case then combustion efficiency would be significantly 
improved by moving the injector further away from the wall. A 
short series of tests has recently been performed using a square 
duct to test this idea. The work has not been documented yet, but 
the preliminary results show a complete lack of burning in the 
larger duct. This tends to indicate that rather than preventing 
combustion, the wall boundary layer may be an important ignition 
aid for the scramjet models used in this study. Some of the 
possible reasons for this are discussed in the section on the 
constant area duct. The reason the boundary layer was ineffective 
as an ignition source may be that the boundary layer compression, 
through the displacement thickness, produced a smaller pressure 
rise in the larger duct. In the small duct the gap between the 
injector and the wall is only 10 mm, allowing time for 
reflections of compression waves within the length of the 
injector. In the large duct this gap is increased to 22.5 mm and 
this would lead to a corresponding drop in pressure increase due 
•to the boundary layer. The model was reassembled with a longer 
injector strut, to check if the extra boundary layer development 
would help with ignition. In this case also no combustion was 
observed. 

Some apparent contradictions are evident in the supposed role the 
boundary layer is playing in the ignition process. In the small 
duct the longer injector strut was found to enhance ignition, 
whilst at the same time wall static pressures in the absence of 
injection were slightly reduced in the constant area section. 
This may perhaps be explained by noting that in supersonic flow 
strong transverse pressure gradients may exist in the flow, and 
that in the immediate vicinity of the injector the pressure may 
not have been reduced with the long strut. The high local 
pressure seen with the internal injector may have been due to the 
wave pattern from the injector leading edge, and might not apply 
on the centre line. Additionally the entrainment into the mixing 
layer of more heated gas may have enhanced ignition. 

The failure of the long strut to induce ignition in the large 
duct is not understood at present. The small duct may have been 
just on the low pressure limit of combustion,' and the reduced 
boundary layer compression in the large duct may have been just 
enough to quench the flame. To check this , the large duct was 
run with a 5 degree compression to increase the pressure and 
reduce the mach number slightly. Combustion was then achieved in 
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the large duct. It would appear that at the pressure levels 
Z'lll / oc !? leva Jle ln the tunnel hypersonic combustion is only 
n?hl Anal ‘ k In . thls condition several effects, which might 

or h Dartial?v° f T 1 SeCOn ? ary > mportance ’ ®ay be able to completely 
°*L_ partial * y prevent combustion. Further testing at higher 

J~r -fJ 1 + 6 ne f essa f y to Properly evaluate the significance 
livers^ v. At thC moment combustion chamber boundary 

k have . a st rong influence on combustion. It 
remains to be seen if this will still be the case when the 
experiments are scaled to larger dimensions and pressures. 

9 ^ he resu lts are presented in the form of specific 
1“^“ ? e against stagnation enthalpy. This shows the hypersonic 

the supersonio l0W temperatura cut out at a higher enthalpy than 
the supersonic, as is to be expected from the lower static 
temperatures associated with higher mach numbers It S!so shols 
that above an enthalpy of about 7 Mj/kg hypersonic combustion 
produces more thrust than supersonic. This is because hypersonic 
flow allows the combustion chamber static temperature to be kept 
in the region of peak specific impulse. The performance at the 
higher enthalpies is well below theoretical maxima, based on 
equilibrium combustion, but it is does demonstrate that at 
hypersonic combustion can under some circumstances be better than 
supersonic. 


Conclusions. 

Tests with a constant area duct show that hypersonic combustion 
is possible with central injection at static intake pressures of 
about 20 kpa. At stagnation enthalpies above 13 Mj/kg, 
corresponding to static intake temperatures of about 1900 K, no 
measurable combustion effects could be observed in comparison to 
the injection of helium at the same mass flow rate. At lower 
enthalpies significant heat release occurs, as evidenced by the 
static pressure rise in the duct. At the higher enthalpies the 
duct static pressure levels agreed well with non equilibrium 
nozzle flow calculations. At the lower enthalpies pressure levels 
agreed with calculations in the upstream sections of the duct, 
but increased further downstream, in some cases doubling the 
theoretical values. This rise is shown to be compatible with the 
presence of turbulent boundary layers in the duct, but should not 
be taken as confirmation of such a condition. 

I n the thrust producing mode values of specific impulse above 
800 sec were measured. Combustion was observed in a diverging 
^° r con< ^i^* ons w ibh static intake temperatures as low as 
400 K, which would require some means of ignition enhancement. 
Viscous heating in the hypersonic boundary layer is proposed as a 
possible mechanism. 

A comparison was made between supersonic (intake mach number 4 to 
4.54) and hypersonic (intake mach number 4.8 to 6.3) combustion 
by means of an intake attachment made of symmetrical inclined 
wedges. Intake static pressure was maintained at the same 
approximate value of 20 kpa. Hypersonic combustion was shown to 
give improved performance at enthalpies above 7 Mj/kg. The 
specific impulse achieved correlated approximately to the intake 
static temperature rather than to mach number. .The exception to 
this is in the low temperature cut off region where the 
hypersonic combustion shows a gradual reduction' in thrust 
compared to the sudden extinction shown in the supersonic case' 
Increased boundary layer heating in the hypersonic flow is 
suggested as a possible cause of this effect. 

32 



9 


The flow conditions currently achieved by the tunnel appear to be 
marginal for hypersonic combustion. Small changes in model 
geometry can cause substantial reduction in the amount of heat 
release in the duct. In this regime the state of the combustion 
chamber boundary layers appears to have a strong influence on 
combustion, but it is not clear yet if this will still apply at 
higher pressures. 
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FIG 1 SCHEMATIC OF SCRAMJET MODEL 
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Fig 2 Constant area duct, helium and tydrogen injection 
stagnation enthalpy 26-5Mj/kg. 
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Fig 3 Constant area duct, helium and hydrogen injection 
stagnation enthalpy &. 39 Mj/kg 
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Fig 10 Comparison of st^er/hyper sonic combustion 
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Fig 11 Computed laminar boundary layers 
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Fig 16 Approximate path of compression waves 
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SHOCK TUNN EL SIMULATIONS . 
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Department of Mechanical Engineering. 
University of Queensland. 

St Lucia. Queensland, 4067. 


ABSTRACT, 


This paper presents the results of 
experiments on a scramjet combustion 
chamber with sidewall injection performed 
in the free piston shock tunnel. T3, at the 
Australian National University. A two 
dimensional model was used with provision 
for parallel and transverse injection. The 
results indicated that combustion was 
strongly influenced bv a region of fuel 
whose temperature was held below its 
ignition temperature by wall cooling 
effects. This leads to a large amount of 
unburned fuel, and produces a significant 
drop ir» specific impulse This is a feature 
of shock tunnel experiments with cold walls 
which would not necessarily applv to a 
flight situation. where aerodynamic 
heating would keep the walls above the 
hydrogen ignition temperature. Maximum 
specific impulse was produced by using a 
combination of parallel and transverse 
injection in a long combustion chamber 
followed by a dual stage expansion. The 
presence of a layer of fuel attached to the 
wall was seen to significantly reduce the 
surface heat transfer rates, indicating the 
potential for weight savings because of 
reduced thermal insulation requirements. 
However, this effect may be dependent on 
the wall quenching layer and might not 
mpply to a flight situation. 


INTRODUCTION, 

Scram jets potentially represent the most 
attractive form of propulsion for high 
altitude and Mach number flight in the 
upper atmosphere, down to the point where 
subsonic combustion becomes more efficient 
at about Mach 4. Ref ill. Due to the nigh 
flight speeds involved, considerable 
aerodynamic heating of the exposed surfaces 
occurs. Active cooling of these surfaces 
will be required in the ducted sections of 
the engine where radiative cooling cannot 
be used Liquid hydrogen is carried on the 
craft as a fuel, and is therefore the most 
suitable coolant. 


up to flight speeds of Msch 10. Above this 
speed sdditionsl hydrogen. or some other 
coolant, must be carried in order to 
maintain a heat balance at reasonable 
surface temperatures. This weight penalty 
would' reduce the useful payload of the 
craft, and therefore engine design must aim 
to reduce heating loads wherever 
possible. 

A particular heat transfer and materials 
problem Is presented bv the construction of 
the fuel injection system itself Fros a 
combustion point of view, an injector strut 
projecting into the flow is desirable, as 
this maximises mixing and fuel heating 
rates by providing two fuel sir interfaces. 
However the total heat input to the strut 
is high because it is heated from two 
surfaces. The strut would be made as thin 
as possible to minimise drag. and this 
would compound the problem of removing the 
heat by means of a coolant. 

Cooling requirements are considerably 
alleviated for the wall injector because 
there is only one heating surface, and 
there is not the same space restriction for 
the coolant passages. At Mach numbers above 
10 hydrogen in excess of propulsion needs 
must be carried, and this will be ejected 
unburned from the craft. If wall injectors 
are used then the excess hydrogen will fcrm 
an unburned layer attached to the wall, and 
this will provide further thermal 
protection. 

The reduced cooling requirements of wall 
injection will be of little use if it does 
not perform well in the thrust producing 
mode. There is now considerable 
experimental data available for injection 
struts, and a basic understanding has been 
obtained about the associated thrust 
producing mechanisms. Ref. (3) has studied 
large scale models up to flight speeds of 
2.2 ka/s in blow down tunnels, and Ref. (4) 
has done chock tunnel studies up to 5.5 
ka/s . 

A program of experimental research into 
wall injection was started in 1985 unaer 
contract to the NASA Langley Research 
Center. The preliminary results of this 
work were presented in Ref. (5), using a 
simple geometry with the thrust surface 
formed behind a single Prandtl Meyer 
expansion fan. Two effects were noticed 
which lead to reduced performance compared 
to central , strut mounted injectors 


Studies have shown . Ref. (2), that the 
hvdrogen required for propulsion is 
adequate to meet the cooling requirements 


♦ s euar Lecturer, Menem c»l Engineer me, Unbar AIAA. 

• Senior fctMWdi Militant. Mechanical Enfineefin*. 

- fcnctreh Militant. Mechanical Emineenne. 

■ Proftiior, Mechanical Eneineerinf, Unbar A1AA. 


Firstly. with the mixing snd combustion 
gone being attached to the wall, net thrust 
due to expanding the combustion products 
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»K ly d * v * i ° p * d « v « r • small section of 
the thrust surface for the short cosbustion 
chaab«r studied. 

Secondly, reduced coebustion mss observed 
end the existence of a quenched lever of 
fuel attached to the wall was postulated as 
an explanation for this effect No 
cosbustion at all could be achieved’ below 
an equivalence ratio of about 1 5 when 
more hydrogen was injected it would burn 
well. but the heat release obtained 
indicated that the hydrogen in the 
immediate vicinity of the wall was still 
unburned. This lead to low values of 
specific impulse. 

In the light of these results a further 

tM “ W * re designed to 

optimise specific impulse within the 

oirti!ii! n * ° f with cold walls and 

partially quenched mixing layers, and to 

gain an understanding of thrust production 
mechanisms in this configuration 

exhau«t n “ ° f t V° ***** expansions in the 
exhaust nozzle. extended combustion 
chambers and a combination of parallel and 
transverse injection it was possible to 
increase the thrust and specific impulse 
developed. However the quenching effect of 
having the combustion zone attached to the 

onlv .iJi ?*? 11 f ° und to Prominent, and 
“* U i*Provements were observed, wall 

°f * U * n *’ Hhic h i* abl* to burn 
spontaneously at room temperature was 

££ thmt th * 1** of combustion 

•!«'"« l.y.r i. not 
expected to be a problem in a flight 

will^fnd « ith .•« rod >'"»*ieally heated 
walls and several means of making shock 

avr^Tihi tMt i nf “ or * representative are 
available. These include the addition of a 
small amount of silane to the fuel 
preheated fuel, and preheated walls. 


EXPERIMENTAL ap paratus 

The experiments were performed in a two 
dimensional model chosen to follow as 
closely as possible the dimensions of a 
previous model with a central injector 
; * • ch#B * tic of the model is shown in 
.* 1 Injection was from behind a 5 mm 

11V * h * ? ,U - The hel « ht of the step 

* ff* n ? lcai t0 the thickness of the 
i, injector strut in order to give 
comparisons between the two 
geometries. The construction of the 

rfrTn X “ Sh °“ n in F1 * 2 Equivalence 
"** , controlled bv means of 
interchangable throat blocks, and injection 


Mach nu«b*r was s«t bv 
reservoir pressure. Throat blocki wire 
available to produce a parallel let Tnd 
Jets inclined at >5 and delr— 1 L 1*1 

f ir 

-"tin 

drilled at *5 degrees to the flow and 
pointing downstream • nc1 

The inlector was followed bv a fatal W 

th * ldn « th 01 which could be 

nozzle of tM th n 25 K * nd 1S ° “ m Th * *xhaust 
nozzle of the combustor consisted of flat 

thrust surfaces, with provision for staged 

expansions and variable divergence angles 

transfer were measured bv 
means of transducers mounted on the intake 
»nd on the thrust surface. The shock 
tunnel was operated in the reflected shock 
mode, with a each 3.5 contoured nozzle 

intake^f conditions «PPropriate to the 
intake of a scram jet combustion chamber 

* PP * r * tu * Procedure is 
d«*cnbed Bore fullv in R«f (*}. 


B I SHANISffg V tL THKUST PWOOUCTinw 

Pr ^ UC !* h * At i r. the 

1 whlch can in turn prndur* 

Th-rl*!« r St ln th * WMion nozzle 
Th.re are two mechanisms bv which thrust is 

, th * burnin * to 

"f* 1 * f th ? ■ od#1 - Compression waves from 
in*fig° 3 ** Ch * ni *** * r * • houn * ch *®sticallv 

riaes^o 1"* h * at r * le *** causes the burned 

111*1 throu-K^; * nd thl * ** nd# com Pression 
waves through the rest of the flow Th«- 

wave. reflect off the model wall, reusing* 
increased local static pressure, 'and f 
corresponding increase in thrust 

Secondly. the fuel Jet creates a r-gion of 

soMd of Ch « nu ^ be ^’ due to the increased 

th^ few Sf-h n Z2 r co “ bu *“c-n. and when 
tnis low Ma^h number zone is sub«»<au-ntlv 

smal ler^Dreas^* n °* Zl * it will expedience a 
r P “ Ur * drop than the high mach 
number freestream or the fuel off flow The 

in**thI*y" 1 *** tCh cr ** tes compression waves 
th. *^*1 freestream and expansion waves In 

ri!ch U TL J#t ;, Uh * n th * co *Pression waves 
reach the wall they create a region of 

thfult*?^ surface pressure. and hence a 
Thi * mechanism of thrust 
Protection is subsequently referred to as 
the expansion interaction effect. it the 
Jet Mach number drops below i <,1 then the 
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effects are reversed, with expansion waves 
in the freestream and compressions in the 
Jet. 


For central injection the low »ach number 
region is in the middle of the duct, and 
the thrust generated by the expansion 
interaction mechanism is realised at some 
distance downstream of the corner. For wall 
in iection the low each number flow is 

attached to the wall. and thrust is 

generated in the region immediately 
downstream of the corner, until the point 
where the pressure differential between the 
burnt gases and the frme*tr*am is 

eliminated. The thickness of the fuel/air 
jet determines the distance over which 
thrust is developed bv the interaction 

mechanism . Wall injection is limited in 
comparison to central injection because the 
mixing layer has only one fuel 
interface. and consequently would not be 
expected to develop so rapidly. Also the 
distance required for the expansion to 
propagate from the corner to the )et 
permits extra mixing with central 
in jec tion 


RESULTS. 


Extra length in the combustion chamber 
before expanding the flow mav be expected 
to allow the mixing layer to spread 
further from the wall across the duct, 
producing more thrust upon expansion of the 
jet. A weight penalty would be associated 
with the longer combustion chamber, and it 
is advantageous to start expanding the flow 
as soon as possible after injection. 
However. expansion of the jet before 
complete combustion is achieved can have a 
quenching effect on the flame, and it also 
slows down the spread oi the jet across 
duct. The optimum configuration of 
combustion chamber and expansion nozzle is 
likelv to be a compromise between these two 
effects . 


It was noted in previous experiments that 
the thrust produced by wall injection was 
not very sensitive to the thrust surface 


jet 

to 

not 


divergence angle when the expansion / 
interaction is absent. This applies 
short combustion chambers which do 
allow for significant propagation of the 
mixing layer from the wall. It 
therefore decided to construct a combustion 
chamber with a 2 stage divergence. 


The first stage includes a thrust surface 
inclined at a degrees to the intake flow. 
This geometry was shown in Kef 5 to produce 
significant thrust due solely to 
compressions from the burning tuel jet. 


Because onlv a smalL pressure drop is 
associated with a a degree expansion . it 
was hoped that the Jet would continue to 
mix and react and spread across the duct in 
the first thrust producing stage. When the 
jet is subsequently expanded by the second 
stage its thickness should be such that 
substantial thrust could be produced by the 
expansion interaction method. 

The advantage over a single expansion of 15 
degrees is twofold. 

Firstly , if the full IS degrees expansion 
is introduced too earlv. no further 
coabustion will take place. By limiting the 
initial expansion to a degrees, thrust aay 
be developed in a region that is still 
supporting combustion. 

Secondly. the final expansion takes place 
at a point where the jet has had time to 
spread away from the wall. and this gives 
potential for increased thrust from the 
expansion interaction effect. 

Short co mbustion chamber. 

Two stage divergen t, 

The initial series of tests was done with 
the model configured as shown in frig 3 with 
a 25 mm length of constant area duct after 
injection before the start ot the degree 
diverging section. 

In Fig a the results of the dual stage 
divergence are compared to previous data 
for a single thrust surface with 15 degrees 
divergence. It is seen that at an enthalpy 
of 8.7 MJ/kg no improvement in performance 
was produced . 

In Fig 5 the pressure against distance 
profiles are shown. It is seen that despite 
significant heat release in the a degree 
section, as shown by the pressure rise 
above fuel off levels, very little net 
thrust is developed on the downstream 
thrust surface. 

At an enthalpy of 4.2 MJ/kg no improvement 
in specific impulse was produced, as can be 
seen from Fig a b. However it can be seen 
from Fig 5 b that the effects of the 
reduced divergence are just beginning to be 
felt in the form of combustion induced 
pressure rises towards the downstream end 
of the first thrust surface. This would 
suggest that somewhat more distance tor 
combustion was required before the start of 
the 15 degree section. Consequently the two 
stage thrust surface was then used with an 
extended combustion chamber , although this 
does to a certain extent defeat the purpose 
of the two stage expansion. which is to 
obtain thrust in all sections where the 
fuel is burning. 
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Lon« co»buitlcn chi»b«r 
4/1S flumi divtr»«ne« . 

At an enthalpy of 8.7 Mj/kf no iaprovenent 
was produced by the extra elxinf length, aa 
eay be seen by coaparing Figs 6a and 4 a for 
specific lapulse. and Figs 7a and 5a for 
P/X dependence. 

At 4.2 MJ/kg a noticeable improvement in 
performance was gained by adding the 

combustion chamber extension. This is shown 
both by the increased value of specific 
impulse . Fig 6b and also by the 

development of net thrust on the second 
thrust surface. Fig 7b. 

It would appear that at the higher enthalpy 
condition .0.7 MJ/kg the temperature 
after the initial expansion of 4 degrees is 
sufficiently high to produce rapid 

combustion . and the heat release is onlv 
limited by the quenched zone attached to 
the wall. The addition of extra combustion 
chamber length does not change this, and no 
increase in the difference between fuel on 
and fuel off was observed. as oiav be seen 
by comparing Figs Sa and 7a. 

However at the lower enthalpy condition 
4.2 HJ/kg it is seen in Fig 5b that 

without the combustion chamber extension 
heat release only occurs 
towards the end of the 4 degree section 
and not much net thrust is produced on 
•i**** r surface. This is thought to be due 
to the longer ignition delay at the lower 
fmpwmtur*. In this case .when a longer 
combustion chamber was used, ignition 
occurred upstream of the first expansion 
and increased thrust was developed on both 
surfaces, as seen in Fig 7b. 

A similar effect has previously been 
observed with central injection, Kef 9 
where increasing combustion chamber length 
ia only beneficial at the lower enthalpies. 
This may also be true for wall injection. 

Trsngverse and paral lel inlection. 

In an attempt to improve the performance of 
the wall injected scramjet a modified 
injector was constructed with transverse 
holes as well as the parallel injection 
throat. It was hoped that the transverse 
momentum of the Jets would carry some of 
the fuel through the boundary layer and 
away from thm qumnching effect of the wall, 
and would also increase the mixing rates . 

A schematic of the injector is shown in Fig 
® * lz consists of a aeries of holes drilled 
•t 45 degrees to the flow and pointing 
downstream. The percentage of transversely 
injected hydrogen was controlled by 


changing the throat of the parallel 
injector. and setting the hydrogen 
reservoir pressure to give the required 
total equivalence ratio. Transverse 
hydrogen mass flow rates equal to 27* and 
79* of that passing through the parallel 
injectors were produced for the o . 9 am and 
the O.l am throats respectively No other 
combinations were used. 

Reduction of the size of the quenched zone 
would be evident in the form of increased 
values of specific impulse, and also 
combustion would be possible at lower 
values of equivalence ratio 

Combustion in a mcramiet using fuel fros a 
room temperature reservoir requires heating 
of the fuel to its 1 ignition temperature tv 
transfer of heat from the free stream, in 
the wall injection case this flow of heat 
is partially offset bv the flow of heat to 
the wall. and is also limited bv or.lv 
having one fuel air interlace. as opposea 
to two for central injectors in tne two 
dimensional configuration the area for heat 
transfer from the flow to the jet is equal 
to the area for transfer from the jet tc 
wall. and this limits the temperature 
which may be achieved in the jet However a 
circular Jet propagating transversely 
across a duct will have a much larger area 
exposed to the flow and mav be expected 
to get hotter . 

Transverse injection 
long combustion chamber. 

9 the results of transverse 
injection with an extended combustion 
chamber are shown. 

Transverse Injection at the 8.7 MJ/kg 

enthalpy condition again gives no 

improvement over any of the other results 
except at low values of equivalence ratio 
Combustion with transverse injection 
occurred at lower equivalence ratios than 
was possible with parallel injection alone 
This is indicated by the two points on Fig 
9a at equivalence ratios of u B3 and l n 
This would suggest that above a certain 
injection pressure the transverse jets 
are to some extent penetrating the boundary 
layer, and burning at lower equivalence 
ratios than is possible for wall injection 
However. the heat release from this 
combustion is not reducing the size of th* 
quenched zone in the parallel injected 
component, and mo no increase in specific 
impulse is observed mt higher equivalence 
ratios. 

Algo mhown on Fig 9s is s dsta point 
obtained by injecting Helium instead of 
hydrogen. This was done for selected 
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conditions, together with the injection of 
hydrogen into nitrogen test gas. to 
separate the effects of combustion from the 
physical presence of a jet of foreign gas 
in the flow. In this case a genuine 

conbustion effect appears to be present, as 
the hvdrogen is producing significantly 
more specific impulse than the helium gas. 

At 4.2 MJ/kg combustion at low equivalence 
ratios is again observed. In addition there 
also appears to be significant improvement 
in performance in the equivalence ratio 
range of 1 to 2.5 , as mav be seen from Fig 
9b. Also shown on this figure are the 

results of 4.2 MJ/kg tests with a 

transverse component equal to 79/1 of the 
parallel injection, the onlv condition for 
which this was done This shows no 

improvement over the 27ft case It had teen 
hoped that a larger proportion of 

transversely injected fuel would lead to 
more combustion. and higher specific 

impulse. This effect would seem to indicate 
that the transverse jets are not 
penetrating verv far into the flow, and are 
still restricted by wall quenching 

Further evidence of this is given in Fig 
10. which compares the wall pressure and 
temperature profiles for a single IS degree 
expansion at an enthalpy of pi MJ/kg tor 
parallel and *7ft transverse injection. No 
difference is apparent between the pressure 
traces . and onlv a slight increase in heat 
transfer with transverse injection was 
observed. Both tests were taken at an 
equivalence ratio ot about 2. where no 
thrust increment was observed with 
transverse injection. A more significant 
result would be at equivalence ratios of 
order one where the transverse injection 
appears to be effective. but no heat 
transfer data was taken at those 
conditions 

It is possible that the transverse jets 
also require a laver of fuel as a thermal 
buffer to insulate them from the wall, and 
if this is not supplied bv the parallel jet 
then more of the transverse component will 
be quenched. Another factor which may be 
significant is the effect of the expansion 
from the trailing edge of the iniectcr The 
strength of this expansion is dependent on 
the amount of fuel injected through the 
parallel throat. Less tuel trow the 

parallel throat would lead to a stronger 
expansion and a region of cooler fuel 
downstream of the injector. 


The presence of the cold model walls 
represents a sink of heat to the flow, 
because in the short duration of the tests 
the wall temperature does not rise 
significantly above ambient. The cooling 
effect of the wall penetrates a significant 
distance into the flow. With central 
injection this does not have a critical 
effect on the development of combustion 

However, when the fuel is injected from the 
wall in shock tunnel testing. there is 
always a region whose temperature wi)l be 
held below the ignition temperature, 
regardless of how much combustion mav take 
place further awav from the wall. The 
hydrogen contained in this low temperature 
region appears to correspond to an 
equivalence ratio of approximately 15. 
because no ignition at all is possible at 
lower equivalence ratios. This represents a 
serious defect in the ability of shock 
tunnel tests to accurately model a real 
flight situation with aerodynamical lv 
heated walls. 

The fuel in the quenched region mav be 
mixed with oxvgen. with onlv trie low 
temperature inhibiting combustion. To 
confirm that this is indeed the case, and 
that the lack of combustion is not. the 
result of some other cause, a test was done 
with the injection of a 20ft silane hydrogen 
mixture . 

The fuel was injected at an equivalence 
ratio of 0.5 into a flow of enthalpy 4,2 
Mj/kg. with hydrogen injection alcne no 
combustion would be expected at this 
condition. The results of this test are 
shown in Fig 11. and it can be seen bv the 
pressure rise above the tuel c-tt levels 
that it burnt well. 

This result is significant in That it 
demonstrates tnat oxvgen is diffusing to 
the fuel laver close to the wall. and that 
it is thermal* effects which are preventing 
it from burning. This gives encouragement 
to efforts which are currently being made 
to design a model which can use heated 
fuel. and possibly even heated walls. It 
also suggests that in the absence cf trie 
above. silane might be used to investigate 
other aspects of combustion which are 
difficult to studv properly in the presence 
of an extensive laver of quenched fuel. 


CONCLUSIONS, 

Significant combustion was achieved in a 
scram jet with wall injection. but the 
presence of cold walls lead to the 
quenching of a significant fraction of the 
fuel, and this produced reduced values of 
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specific impulse This effect was partially 
offset bv the use of transverse injection 
and staged expansions. The successful 
combustion of a silane hydrogen mixture at 
an equivalence ratio of 0.5 was taken as 
confirmation that the reduced combustion 
was due to thermal effects, and was not the 
result of some other limitations of the 
facility . 

To make shock tunnel testing with wall 
injection more representative of flight 
conditions where hot walls would not be 
expected to quench the fuel . the use of 
preheated fuel. preheated walls or the 
injection of small amounts of silane are 
being investigated. 


ACKNOWLEDGMENTS 

This work was pertormed under a grant t rom 
the NASA Langley Research Center. 
Hypersonic® Propulsion Branch 

The experimental work was performed in the 
shock tunnel ot the Physics department of 
the Australian National University, and the 
assistance of manv members ot the statt tor 
all the back up support required for the 
operation of a major test facility is 
greatlv appreciated. 


T ABLE OF TEST COftPJTiPN S 

H M T P 

MJ/kg K kPa 


0.70 

3.50 

2500 

160 

A. 20 

3. SO 

1100 

1 60 



1. Robert A. Jones. Paul w. Huber. "Towards 
Scram Jet Aircraft". 

Technology Report. AIAA Journal of 
Astronautics and Aeronautics. 

Feb. 1978. 

2. Griffin Y. Anderson. "An Examination of 
In jec tor /Combustor 

Design effects on Scram jet Performance" 
Presented at the 2nd International 
Symposium on Air Breathing Engines. 
Sheffield, England. March 1974. 

3. Robert w. Guv. Earnest A. Macklev. 
"Initial Wind Tunnel Tests at 

Mach l and 7 of a Hydrogen burning- 
Airframe integrated Scramjet" 

Presented at the 4th international 
Symposium on Air Breathing Engines. 

Lake Buena Vista. Florida. April 1-6. 1979. 



IXTJlf FLO# 



conmim mmu anisic* imincriot imcis 

conntssim m atmim coamssim m annum 


INTAKE 



*mi lira 


4. R . J . Stalker . R .G Morgan 

"Supersonic Combustion with a Short Thrust 
Nozzle". 

Journal of Combustion and Flame. Vol. 57. 
No. 1 . 

July 1984. pp 55-70. 

5. R . G . Morgan . A . Paul 1 , 

N. A Morris .R. J . Stalker 

"Hydrogen Scramjet With Sidewall 

Injection” 

Paper presented to the 2nd l.E.Aust. 
National Symposium on Space 

Engineer mg . Sydney Mav 1986. 


Fif ». fe*l tumd teemion Cttfffl UMKSI* 



Fit 31. Came e*vmiei/«ll M intartctiec 


51 


ORIGINAL PAGE IS 
OF POOR QUALITY 



specific 
impulse (seel 


stagnation enthalpy 97 Mjlkg 


500 

— i i — i — ' » — » — 


400 ■ 

a 2 stage 4/75* 


o o 1 stage IS * 

* 

300 ■ 

□ 



O Q 


200 

l 

0 

100 

o 0 Q o 

* 

c 


15 


equivalence ratio 

dual stage 

Fig 4 A. Comparison of sin^e and 
expansions. Parallel wal l injection 


short combustion th»iber. 


specific 


impulse (s*cl stagnation enthalpy U2 Mjlkg 

500 

» i * » ■ r * 

* » 

400 

a 2 stage 4/15* 
o 1 stage 15 9 

- 

300 

m 

a 

200 

: V • 

* 

100 

o o 

0 





' ^ 1 ' 5 

15 


♦ 

0 


0 02 U 


001 


intake pressure 
‘ *■ stagnation 
pressure 


002 


equivalence ratio 
Fig kd . Comparison of sin^e and dual stag* 
expansions. Parallel wal l irijeitioa 
short combustion chamber. 

intake pressure legend for eqji valence ratio 
* stagnation ° P * 2-1® 


n 1 rr - 

— 1 1 

15° 

J 4° ' 

- i 


|09 0 9 Qjng 1 


‘ p °° 0 0 0 ! 
f °oO°P 

l 

its, i 

1 l L_ 

a— 1 1 


legend for equivalence ratio 
n 0 =162 
o 0 

♦ helium ir 


0 100 200 300 400 

distance from intake (mm) 
stagnahon enthalpy 67 Mjlkg 


500 600 


F»« SA. Dual stags axpansion, short duct, 
Parallax wall injaction. 



ooi r 


0 100 2 00 300 400 

distance from intake (mm) 
stagnahon enthalpy 4 -2 Mj/kg 

Fie SB. Dual stags axpansion, short duct. 
Parallel «all Injaction. 


specific 
impulse (sec) 
500 



0 12 3 

equivalence ratio 

Fig 6 A. Stagnation ^enthalpy 8-7 Mj/kg 
legend 


long duct 
o short duct 


specific 
impulse (sec) 
500 


stcgn,hm jmHV 



0 12 3 

equivalence ratio 

Fig 6B.0ual stage expansion, long and short 

combustion chambers. Parallel injection. 


52 









ORIGINAL PAGE IS 
Of POOR QUALITY 


static pressure o equivalence rati o 2 kU 
stagnation o 154 

pressure * fuel off 



Fig7B stagnation enthalpy 8 7 Mj/kg 


, o equivalence ratio 182 

static/ ♦ V18 

stagnation pressire o fuel off 



stagnation enthalpy U-2 Mj/kg 
Fig 7A Dual stage expansion, long combustion 
chamber, parallel wall injection. 


trinsmt 



specific 


impu 

1000 

tse (sec) 


""I 1 * “ 

i 

800 

• 

- 

600 

$* * 

- 

4 00 

• o • . 

o o 


200 

• • 

a ■ 

i ■ * 

■ 

« 


i i 1 1 

0 1 2 3 4 S 


equivalence ratio 
B stagnation enthalpy 4 2 Mj/kg 

O 2 Stiff i#H 4 MCt,??/atm*V 8 r$t Jfl/fCtjOA. 

• 2 stiff nNMion. 

O 2 stiff fiPWfiM' Jo»f toct, Hnllil j»;rctio*. 

• 2 stif # nwiuott, *A*rt *»ct. H nlifl miKtio*. 

• J Stiff flPIMlO*. S**ft IfCt. JAJKtJO*. 

m j stiff Mellon, ftort *tct'27%trM»%*nt, hthvt 


specific 
impulse (sec) 

f/t/l - 


1 1 1 * 

4 00 

m m 

■ 

300 

• 

♦ «. *o 

200 

V • 

100 

♦ * ■ . 

■ 

u 1 1 


0 1 2 3 4 5 

equivalence ratio 
A. stagnation enthalpy 8 7 Mj/kg 


Fig 9 Various configurations of 
wall injection. 


53 







static 

pressure kPa 



distance from intake (mm) 


A. Pressure profiles 

a equivalence ratio =2 06 
wa ^ ■ 27%transverse injection 

temperature(K) * equivalence ratio = J W 
ca parallel injection 



o 10 0 200 3 00 400 500 600 


distance from intake (mm) 

B. Temperature profiles (transient) 

Fig 10. Comparison of parallel and transverse 


wall injection, long duct, IS 8 expansion 


static/ o fuel off 



0 100 200 300 400 500 600 700 

distance from intake (mm) 


Fig 11. Wall injection of silane/hydrogen 
mixture Constant area duct 



54 






NUMERICAL MODELLING OF SIDEWALL INJXCTZD SCRAMJETS 


C . BraacianinI t 


R.G. Morgan 


55 




NUMERICAL MODELLING OF SIDEWALL INJECTED SCRAMJET 


Side wall injection is an attractive alternative to central injection 

in a scramjet engine, since the hydrogen layer can shield the walls of the 

scramjet engine from the high temperatures found in the free stream. 

Previous experiments using a model scramjet in the high enthalpy flow 
created by shock tunnels, however, have indicated that the performance of 
the wall Injected scramjet is significantly reduced when compared with the 
performance of the central injected scramjet (Ref. 1). A numerical study of 
the scramjet's flowfield was undertaken to see if the experimental results 
could be reproduced; and if they could, to see what other details of the 

flow might be learned from the numerical results. 

All scramjet experiments in the shock tunnel, to date, have been 

performed using fuel with a total temperature equal to the local room 
temperature. In future experiments it has been proposed to heat the 
hydrogen fuel so that the temperatures will better simulate those which 
occur in a real flight vehicle. The numerical program was also used here as 
a predictive tool, so that the likely changes in performance due to fuel 
heating may be ascertained. 

DESCRIPTION OF COMPUTER PROGRAM 

A two-dimensional, parabolic computer program was used to compute the 
flow field inside the scramjet model. The computer program (known as CHARMS 
at the University of Queensland) is based on the CHARNAL computer program 
which is described in Ref. 2. The program is specifically designed to 
calculate the turbulent mixing and chemical reactions which occur when 
hydrogen is injected parallel to a main stream of air. The program solves 
the time-averaged parabolic partial differential equations for the 

transport of momentum, energy, and species. The solution to the 

differential equations is accomplished by the finite difference method of 
Patanker and Spalding (3). 

A finite rate chemistry scheme is used in the program to calculate the 
chemical reaction rates. This scheme is described in more detail in Ref. 4. 
This reference describes the use of two available reaction systems. The 

first involves the use of 7 species and 8 reactions (treating N^ as inert), 

and the second Involves 12 species and 25 reactions. In addition the 
chemical reaction rates may be reduced to allow for the effects of 
"unmixedness” often found in turbulent reacting flows. For the cases 
studied in Ref. 4 it was found that the 25 reaction system was superior in 
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predicting Ignition, however once Ignition had occurred the 8 reaction 
system was as good as the 26 reaction system. For all the cases studied In 
this report the 8 reaction system has been used, and the effects of 
unmlxedness have not been included. 

The turbulent viscosity In the program is calculated by the well known 
k~c model of turbulence (Ref. 6). In Its most common form (known as the kei 
model) the viscosity is found from the formula: 

A = C u <1> 

H c 

where the quantities k and t are found by solving a pair of transport 

equations simultaneously with the equations governing the mean flow, and C 

. » 

is a constant. 

The ke2 model is an extended version of the k-e model which contains a 
weak shear flow correction (Ref. 6) in the form of 

C H = 6<p7Z> (2) 

Here, P/e represents the average value of P/e across the layer. This 
correction greatly improves the k-e models's ability to predict flows where 
the production P and dissipation e of the turbulence are not in balance. 

The k-e turbulence model was developed for incompressible flows, and 
its application to compressible shear layers has shown that it tends to 
overestimate the mixing rate in high Mach number flows. To help overcome 
this problem an empirical compressibility correction has been developed 
(Ref. 6) which is applied to the turbulent viscosity calculated from the 
k-e model. This correction factor is applied whenever the Mach number is 
above one. The correction factor is evaluated as follows: 

K(M t ) = 0.25 + 0.75 / ( 1.0 + exp( 24.73( M^ - 0.2 )) (3) 

where M^ is k 1/2 divided by the local speed of sound. 

The compressibility corrected version of the ke2 model has been used 
to obtain all the results presented In this report. The turbulence 
constants used In the calculations are as follows: 
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C = 0.09 g(P/e) , C = 1.44 , C p _ * 1.92 . x * 0.435 

Cl 

a t k = 1.00 ,-a t c « 1.30 . other 0.9 . a |f# = 0.7 


CHARNAL makes use of wall functions to relate the fluxes through the 
walls with the values of the dependant variables at the near-wall nodes. 
The use of wall functions Is required because the k-c model employed is not 
valid in the low Reynolds number region near the walls. The wall functions 
used in this report assume uniform shear stress prevails in the region near 
the walls, an assumption which Is not entirely valid If large pressure 

gradients exist. Wall functions which do consider the effects of pressure 

gradient are available in CHARNAL, however they have not been used to 

obtain the results presented here due to instability problems. 

Pressure gradients In the lateral direction are calculated in the 

program by using the SIMPLE (Semi-Implicit Method for Pressure-Linked 
Equations) algorithm. The application of this procedure for use in the 
CHARNAL program is described In Ref. 7. 

TEST CONDITIONS 

Pree Stream 

Two nominal , free stream stagnation enthalpy conditions were used, 
namely 4.2 MJ/kg and 8.7 MJ/kg, and one Mach number, M = 3.5 . The values 
of velocity, temperature, and dissociation at the exit of the shock tunnel 
nozzle (and thus at the entrance to the scramjet) were estimated by using 
NENZF (Ref. 8) , which takes into consideration the non-equilibrium effects 
occurring as the flow in the nozzle expands. It should be noted that only 
nominal free stream conditions were used in making the calculations. The 
results are shown in Table 1. 

TABLE 1 FREE STREAM CONDITIONS AT SCRAMJET INLET 


H 

• 

(MJ/kg) 

T 

(K) 

PIN 

(kPa) 

U 

(m/s) 

a 

<%> 

M 

8.7 

2520 

160 

3250 

10.3 

3.34 

4.2 

1160 

158 

2380 

0.18 

3.59 
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Hydrogen 


Hydrogen reservoir (injection) pressures were taken directly from the 
experimentally recorded reservoir pressures where available (room 
temperature fuel). When the experimental hydrogen reservoir pressures were 
not available (heated fuel runs), the same reservoir pressures as the 
equivalent room temperature run was used. 

INITIAL CONDITIONS FOR COMPUTER PROGRAM 


Hydrogen static pressures at the exit of the injector were generally 
not matched with the free stream static pressures during experiments. In 
reality this would result in a strong wave/shock structure immediately 
downstream of the Injector. Current attempts to model this structure with 
CHARNAL have shown stability problems which are yet to be overcome. As a 
result a much more approximate treatment of the flow characteristics near 
the injector have been taken in this report. First the hydrogen and air 
mass flow rates, total temperature, and stagnation pressures are 
calculated. The two gases are then allowed to expand in a l-D isentropic 
fashion, keeping the static pressures the same in both streams, until the 
combined thickness of the two streams match the height of the duct after 
the injector. The velocities, and temperatures of the two streams, and the 
common static pressure, are then used as the initial conditions for the 
program. A step change in the velocity and temperature is assumed at the 
mixing layer. If any recirculating regions are present near the injector, 
their effects have been ignored in the calculations. This is necessary 

since CHARNAL. being a parabolic program, cannot predict regions with 
recirculation. 

A fuel preheating rig is nearing completion at the Mechanical 
Engineering Department. This rig will quickly heat the hydrogen fuel by 
burning some of the hydrogen fuel with oxygen before injection. Initial 
estimates show that the total temperature of the fuel after this heating 
may be as high as 1800 K. with static temperatures of the order 1000 K. For 
this reason hydrogen total temperatures of 1800 K have been used in the 
heated fuel computations. The experimental rig will also result in some 
water vapour being injected along with the hydrogen. CHARNAL currently 
cannot handle the initial presence of water with the injected hydrogen. As 

a result the effects of the initial water vapour have also been neglected 
in this report. 

initial values of turbulent kinetic energy and dissipation length 
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scale are also required before computations can begin. Since these values 
are not measured an estimate Is required. The values used In this report 
are: 

k 

Hydrogen — * 0.005 
U z 

Air = 0.001 

U 

e 

— = 0.03 
DJ 

A constant temperature is used as the boundary condition for the 
energy equation along the top and lower walls. The temperature of the lower 
wall was fixed at 291 K for the calculations. 

RESULTS 

(A) CONSTANT AREA DUCT RESULTS 

H * 4.2 MJ/kg . <P - 1.43 . 

s 

Room Temperature Fuel 

Shock tunnel experiments in Ref. 1 .using the scramjet model with wall 
Injection, had shown that very little pressure rise occurred in the 
scranvjet duct when the equivalence ratios were below approximately 1.6 . In 
an attempt to explain this phenomena It was proposed that a region of 
hydrogen fuel, which was injected along the lower wall of the model, was 
being quenched by the cold walls, and thus preventing combustion from 
occurring. Since the particular test case under study here was also at a 
low equivalence ratio, it presented a good case to test this theory against 
the results from numerical computations. 

Fig. 1 displays a plot of the computed pressure (normalized against 
Initial pressure) versus distance from the injector. Also plotted on the 
same figure are two sets of experimental data taken from Ref. 9. The first 
set are the results when a cylindrical nozzle was used in the Injector, and 
the second are the results when a source-flow nozzle was used. (The shape 
of these nozzles is explained In more detail In Ref. 9). Both sets of 
experimental data have been normalized against the experimentally recorded 
pressure in the inlet to the scramjet (upstream of the Injector). 
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Fig. 1 displays a reasonable amount of agreement between theory and 
experiment when the^ inaccuracies of the experimental pressures, and the 
large assumptions required to obtain the initial conditions for the 
computations are considered. The normalized pressure does not rise above 
1.2S until approximately 27 cm downstream of the injector. 

Fig. 2 shows the experimentally recorded heat transfer results along 
the lower wall of the model (these results have been taken from Ref. 9 and 
have been amended by a factor of 2 ), along with the numerically computed 
results. Again, considering the large scatter in the experiments, the 
agreement is very good. 

Since there is a reasonable amount of agreement between the 
computations and all the experimental results currently available at this 
particular test condition, it is now interesting to turn to the 
computations to see what other information may be obtained. 

Fig. 3 displays the mass fraction of water versus distance from the 
lower wall at 3 different distances downstream of the injector. It can be 
seen that very little reaction has occurred at a distance of 10 cm. 

Temperature contours are shown in Fig. 4 .with distances downstream of 
the injector along the lower side, and the 2.5 cm duct height along the 
vertical side. Very little temperature rise is seen until approximately 
12 cm downstream. It would appear that the ignition is quite slow at this 
lower enthalpy, delayed by both the low temperatures and low oxygen 
dissociation levels. A 25 reaction chemistry system, if it had been used, 
may have produced a more accurate estimate of the ignition delay length 
than the 8 reaction system used here. 

The temperature contours in Fig. 4 also indicate quite a large cold 
region of hydrogen near the lower wall which could possibly prevent 
combustion. However, it is also important to consider the local mixture 
equivalence ratios. Self ignition may be expected to occur in region where 
the mixture equivalence ratio is approximately 0.2 .and the temperatures 
are above 800 K (Ref. 10). Fig. 6 displays local equivalence ratio 
(ignoring 0 2 and H 2 in the form of water) contours at levels of 0.2, 1.0 
and 1.8 .These contours give some Indication of where the flame front is 
likely to occur. The <P = 0.2 contour can be seen to be very close to the 
T = 750 K contour for some distance downstream of the injector, helping to 
explain the long ignition delay. However, after ignition, the flame front 
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continues to penetrate into the hot main free stream air. The cold 

hydrogen, located near the vails. Is quite a distance from the flame front 
and Its effects on the combustion seem small. The combustion efficiency of 
the mixed hydrogen <t) rr ) after 40 cm was 71% .while only 7% of the total 
hydrogen (0 TF ) had burnt. It would appear that the amount of combustion 
which occurs after Ignition is mainly limited by the rate at which the 
hydrogen Is penetrating Into the free stream. 

Hot Fuel 

The simulation at 4.2 MJ/kg was repeated using hot hydrogen fuel. 

Since the total amount of fuel Injected was identical, but the density was 
lower, the initial thickness of the hydrogen Jet was thicker than for the 
room temperature case, and the initial static pressures where slightly 

higher. 

Normalized pressures versus distance is shown in Fig. 6 ,and indeed 

indicates a larger pressure rise in the duct than was calculated for the 

room temperature hydrogen results. Fig. 7 shows that, unlike the room 
temperature fuel case, ignition has already occurred well before x = 10 cm. 
Water mass fraction contours .in Fig. 8 .show some water forming almost 
immediately after injection. Temperature contours, in Fig. 9 .show that the 
hot regions of the flow extend almost the entire way to the lower wall, 

indicating virtually no effect due to wall quenching. Fig. 10 shows that 

the flame front has penetrated very little into the main stream flow. 

Combustion efficiencies (Fig. 11) show that the combustion efficiency of 
the mixed hydrogen (n RR ) is 92%. while only 15% of the total amount of 
hydrogen available for combustion (rj _„> has been converted to water. 
Again the total amount of hydrogen which has combusted seems to be limited 
mainly by the mixing rate. 

8.7 MJ/kg . <t> = 2.64 

Room Temperature Fuel 

Since the equivalence ratio is well above 1.5 this is a condition 
where a reasonably large pressure rise would be expected. Fig. 12 displays 
the normalized pressure versus distance from the injector. The experimental 
values have again been taken from Ref. 9. The computed pressures seem to 
rise periodically under the influence of some pressure wave travelling down 
the duct. The overall pressure rise predicted is in reasonable agreement 
with experiment, but the location of the • pressure rises is in some 
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question. Again, this Is not surprising considering the Initial condition 
assumptions, and the scatter In the experimental results. The source of the 
pressure wave is evident from the pressure contour diagram shown in 
Fig. 13. Here a strong pressure wave Is seen to emanate from the initial 
mixing region between the hydrogen and the free stream air. This pressure 
wave 16 due mainly to combustion. This was shown by running the program 
again with identical conditions .but without chemical reactions. The 
pressures along the lower wall when reactions where not used is also shown 
on Fig. 12. 

Water mass fraction contours, shown In Fig. 14, show that with 
chemical reactions. Ignition has been almost immediate. No advantage would 
probably have been gained here if the 25 reaction system had been used. 

Fig. 15 displays the predicted heat transfer results, along with 
experiments. The experimental results show a lower heat transfer rate (and 
thus better Insulation for the wall) than predicted by the computer 
results. This may suggest that a lower dissipation length scale should have 
been used in the computations. The heat transfer rates were also found to 
be quite sensitive to the initial conditions chosen. For example, reducing 
the velocity of the free stream air in a small region near the injector, to 
the same velocity as used for the hydrogen (while keeping the total 
temperatures the same), resulted In the heat transfer results shown in 
Fig. 16. Results here are in much closer agreement with experiment. 

Temperature contours shown in Fig. 17 show a large cold region near 
the wall, but this appears to be mainly due to the large amount of cold 
hydrogen which has been injected rather than due to the quenching effects 
of the walls. From Figures 17 and 18 the combustion region can be seen to 
easily extend Into regions where the temperatures are well above those 
required for combustion. The combustion region, however, only extends to 
approximately halfway across the scramjet duct. Fig. 19 shows that 90* of 
the mixed hydrogen available for reaction has reacted (n RB >. however, only 

23* of the total hydrogen available for reaction has been converted to 
water ( t i ). 

STOICH 

Hot Puel 

Figure 20 shows normalized pressure versus distance when heated 
hydrogen is used. Again the initial pressures and Jet thickness are larger 
here than those used in the room temperature case. The total pressure rise 
recorded is 1.36, compared with 1.68 for the room temperature hydrogen. The 
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pressure wave emanating from the Initial mixing region seems to be of a 
much lower strength, than was seen in the room temperature result. 

Temperature contours can be seen in Fig. 21. Figure 22 shows that the 
flame front is spreading into the oxygen rich region, but that the 
spreading rate is quite slow. Combustion efficiencies in Figure 23 show 
that 87% of the mixed hydrogen available for reaction has burnt <n BR ). 
while 18% of the total hydrogen available for reaction has been completely 

converted to water (n sT0ICH )- 

(B) 16* DIVERGENCE 

H =4.2 MJ/kg • <J> = 1-43 
s 

The conditions at 4.2 MJ/kg were repeated except that a 15* diverging 
thrust surface was located on the lower wall, 20 cm downstream from the 
Injector. The thrust surface was extended for 20 cm in the axial direction. 
Both room temperature and hot hydrogen runs were performed. The results for 
room temperature fuel are displayed in Fig. 24, while those for heated fuel 
are shown in Fig. 25. Also shown on these diagrams are the results obtained 
when no chemical reactions where used, and the result when no fuel at all 
was Injected. The pressures have been normalized against the nominal 
pressure in the inlet to the scramjet, before the Injector. This Inlet 
pressure is common to all the computational results presented at this 

stagnation enthalpy. 

By comparing the results with and without reactions for the room 
temperature hydrogen it can easily be seen that ignition is delayed to some 
8 cm downstream of the injector. Pressures then begin to slowly rise above 
the no-reaction case. Approximately 2 cm downstream of the expansion corner 
a slight hump can be seen in the pressure on the thrust surface. This may 
be due to the compression waves produced when an expansion fan from the 
corner interacts with the Mach number gradient in the flow (Ref. 11). 

In comparison the hot hydrogen fuel results (Fig. 26) show almost 
almost Immediately ignition. The pressures obtained when reactions are 
included quickly rise above those obtained without reactions. Again a 
slight hump in pressure is seen in the thrust surface Just downstream of 
the expansion corner. 

A summary of the thrusts obtained are shown in Table 2, the specific 
impulses in Table 3, and the combustion and mixing efficiencies in Table 4. 
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rom Table 3 It can be seen that the specific Impulse for the hot hydrogen 
case is significantly larger than those obtained with the room temperature 
However, note that the gain in specific impulse due to chemical reaction Is 
of the same approximate size In both cases. It would appear that the main 
gain In specific impulse is due to the larger Initial static pressure. 

Table 4 indicates that only 64* of the mixed hydrogen available for 
reaction had reacted by the end of the duct when room temperature hydrogen 
was used, while 93* of the hydrogen had reacted when heated fuel was used. 
It seems that although significantly more hydrogen had reacted when heated 
fuel was used, the gain in pressure was not proportional. 


TABLE 2 COMPUTED THRUST 
H s = 4.2 MJ/kg. * = 1.43. 16- DIVERGENCE x = 40 cm 


THRUST (N) 

ROOM TEMP 

HOT 

WITH REACTIONS 

2461 

3816 

NO REACTIONS 

1849 

3141 

NO FUEL 

1463 

1463 


TABLE 3 COMPUTED SPECIFIC IMPULSE 
H s = 4 2 MJ/k 8 ♦ = 1-43 16 * DIVERGENCE x = 40 cm 


I ( sec ) 

ROOM TEMP 

HOT 

WITH REACTIONS 

267 

400 

NO REACTIONS 

194 

330 


i 
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TABLE 4 COMPUTED EFFICIENCIES 
H s - 4.2 MJ/kg ♦ = 1.43 16* DIVERGENCE x = 40 cm 


T) (%) 

ROOM TEMP 

HOT 

^RR 

54 

93 

^TF 

4 

7 

^STOICH 

6 

11 

Vx 

11 

12 


DISCUSSION 

Numerical results for the test conditions here seem to indicate that 

the finite rate reactions are fast, and that combustion is mainly limited 
by the mixing rate. It is, however, significant to note that at least one 
set of experiments reported in References 1 and 9 would seem to indicate 
that this is not the case. Namely, when the sidewall injected model 
scramjet was configured as a constant area duct, with free stream test 

conditions H^ 4.2 MJ/kg. M * 3.6, and ♦ < 1.3 .very little pressure rise 

was observed in the duct. However, when the same model and test conditions 
were used, but a 20% sllaneiSiH^l/hydrogen mixture (based on mole 
fractions) was injected at $ = 0.69, a significant pressure rise was 

recorded. The fact that the silane/hydrogen mixture produces such a large 
pressure rise would seem to Indicate that there was significant mixing and 
combustion with the free stream. This large pressure rl6e may. in fact, be 
due to the smaller Ignition times, or some other gas dynamic effect which 
has not yet been considered. Experiments to check for gas dynamic effects 
in the central injection configuration have been performed by injecting 
argon gas to simulate the molecular weight of silane, but avoiding the 
chemical reactions. These argon experiments did not produce the large 
pressure rise that was seen when silane was used. This demonstrated that, 
at least in the central injection configuration, the pressure rises were a 
genuine combustion effect. Computer simulations of the silane experiments 
to date have only been one-dimensional, so no simulations which Include 
mixing rates have yet been performed. If these were available it would 
provide a useful comparison for the computer simulations presented here. 

It is also worthwhile mentioning certain points about the hydrogen/alr 
simulations performed in this report which have not yet been considered. 
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Firstly, the effects of unmlxedness on the reaction rates have not been 
determined. If unmlxedness vas taken Into consideration It may Indicate 
that there are certain slow reactions which are significant In delaying the 
combustion rate. Secondly, a single value of In the turbulence model was 
taken across the flow (although varied with axial distance). It may well 
be that a high- turbulence production rate in the boundary layers was 
adversely affecting the calculation of a suitable value of for use In 
the mixing layer. If was allowed to vary across the flow this may tesult 
In a faster mixing rate In the free stream, and a slower diffusion rate In 
the boundary layers. Thirdly. It Is obvious that the method of calculating 
the heat transfer rate to the walls must be reasonably accurate In this 
work if the hypothesis of a wall quenched hydrogen layer is to tested. A 
low Reynolds number version of the k-c model, which does away with the need 
for wall functions, should be able to predict the heat transfer rates with 
greater accuracy. However, such turbulence models are computationally more 
expensive and would require significant alterations to CHARNAL. In any case 
the predicted heat transfer rates here are in reasonable agreement with 
experiment, at least at the 4.2 MJ/kg condition. Fourth, only nominal free 
stream conditions have been used in this report. It would be preferable if 
NENZF, along with the experimentally recorded pressures at the inlet to the 
scramjet, were used to calculate the free stream conditions for each 
individual experimental test case. Finally, the assumption of equal static 
pressures in the hydrogen and air jets at the initial station is obviously 
incorrect. If the current instability problems can be overcome, these 
simulations should be repeated Including the effects of the unmatched free 
stream and hydrogen pressures at the injector, so that the significance of 
the wave/shock interactions can be determined. 

CONCLUSIONS 


The numerical results for pressure and heat transfer In a constant 
area duct, using room temperature fuel, indicate a fair agreement with 
experiment. The numerical results show that at the low enthalpy (4.2 MJ/kg) 
and low equivalence ratio ( <t> = 1.43) case there is a significant ignition 
delay. In comparison the 8.7 MJ/kg and high equivalence ratio ($ = 2.64) 
case showed almost immediate ignition. For both enthalpy conditions the 
flame front spread outwards into the hot. oxygen rich, free stream. The 
cool layer of hydrogen near the walls was due to the large amount of 
hydrigen which had not reacted, or mixed with the free stream, rather than 
due to cold model walls. The total amount of mixed hydrogen which had 
combusted to form water was found to be quite high, while the percentage of 
the total amount of injected hydrogen which had reacted was small. This 
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indicated that the main limiting factor, after ignition, seemed to be the 
rate at which the hydrogen mixed with the free stream, rather than the 
finite rate chemistry. 


Heated hydrogen fuel seemed to have its greatest benefit at the lower 
enthalpies where ignition lengths were greatly reduced, and the total 
amount of reacted hydrogen was significantly increased. However, the 
pressure rise and thrust due to tombustion effects were not significantly 
different. The predicted total thrust and specific impulse of the heated 
fuel scramjet was much improved over the room temperature fuel case. The 
major part of this improvement seems to be due to the larger static 
pressures immediately after the injector, which are due to a gas dynamic 
effect, rather than due to combustion. 


It is well recognized that the k-e turbulence model does have certain 
limitations, however the CHARNAL computer code has performed reasonably 
well considering no adjustment of the turbulence constants were required to 
produce any of the results presented here. Work is currently underway In 
the Department to carry out flow visualizations in the scrarajet model, and 
to measure species profiles. This information will greatly assist in 
providing initial conditions for the CHARNAL program, and should enable a 
fine tuning the turbulence constants, and a more accurate selection of the 
dissipation length scale. 
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SYMBOLS 


V C er c £2 

DJ 

H s 

I 

k 


t 

t 

M 

“r 

P 

PHI 

PIN 

PO 


Constant coefficients appearing In turbulence model 

Injector step height 

Stagnation Enthalpy 

Specific Impulse 

Kinetic energy of turbulence 

Dissipation Length Scale C k 3/2 / c 

Mach number 

Mach number of turbulence 
Production rate of turbulence energy 
Equivalence Ratio ($) 

Nominal static pressure In inlet to scramjet 
Computations: Initial static pressure used in computations 
Experiments: Experimentally recorded pressure in inlet to 
Scramjet 
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q Heat transfer rate 

T Static Temperature 

.T Hydrogen Total Temperature 

TH Thrust 

U Velocity 

x Axial distance downstream of Injector 

y Distance normal to lower wall 

GREEK SYMBOLS 


a 

e 


n 

n 

'mix 


n 


RR 


n 


STOICH 


'TF 


e 


D 


X 


P 


a 

t 

o 

1 


♦ 


Percentage mass of oxygen that has dissociated 
Turbulence energy dissipation rate 
Efficiency 
Mixing Efficiency. 

Defined as the amount of reacted H 2 (H ? in the form of 
water) if all mixed hydrogen and oxygen reacted completely, 
divided by the same quantity If mixing had been 
complete 

Reaction Rate Combustion Efficiency. 

Defined as reacted H 2 .divided by the amount of reacted H 2 if 
the hydrogen and oxygen which are mixed reacted completely 
Stoichiometric Combustion Efficiency. 

Defined as reacted H ? .divided by amount of reacted H 2 if 
mixing was complete and the hydrogen and oxygen reacted 
completely 

Total Fuel Combustion Efficiency. 

Defined as reacted H 2 .divided by the total amount of H 2 

Thrust surface (lower wall) divergence angle 

von Karman’s constant 

Turbulent viscosity 

Density 

Turbulent Prandtl/Schmidt No. 

Laminar Prandtl/Schmidt No. 

Equivalence Ratio 
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"*• 1 Axial variation of praaaura In a conatant araa duct 
♦ - 1.43 . H s - 4.2 MJ/ltJ . II » 3.6 . T, 300 K 



««• 2 Axial variation of h.at tr.n.f.r . „ nttmnt „„ duct 

♦ - 1.43 . H s - 4.2 MJ/kg . 11 « 8.5 . T . 300 K 

71 





H20 MASS FRACTION 



Pig. 3 Profiles of water mass fractions at 3 axial locations in a 
constant area duct 

(J, = 1.43 . H s = 4.2 MJ/kg , If * 3.6 , T t = 300 K 
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Pig. 4 Temperature contours (K) In a constant area duct 
♦ = 1.43, H s = 4.2 MJ/kg . M * 3.6 . T^= 300 K 
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Pig. 6 Local equivalence ratio contours In a constant area duct 
* = 1.43, H g = 4.2 MJ/kg . 14 ■ 3.6 . T t = 300 K 



Fig. 6 AxUl variation of pressure In a constant area duct 
♦ » 1.43. H s « 4.2 MJ/kg . II « 3.6 . 1800 K 
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Pig. 8 Water mass fraction contours In a constant area duct 
♦ = 1.43. H g * 4.2 MJ/kg . M = 3.6 . 1800 K 
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Fig. 9 Temperature contours (K) In * constant area duct 
♦ = 1.43, H g = 4.2 MJ/kg . M = 3.6 , T^= 1800 K 
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Pig. 10 Local equivalence ratio conto 
* = 1.43. H s = 4.2 MJ/kg . M = 3.6 . T 




*' 11 Axial variation of combustion and mixing efficiencies 
* * 1.43. H s = 4.2 MJ/kg . M » 3.6 . T t = 1800 K 
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Fig. 16 Axial variation of heat transfer rate in a constant area duct 
<t> = 2.64 . H s = 8.7 MJ/kg . li = 3.6 . T t = 300 K 
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Fig. 16 Axial variation of heat transfer rate In a constant area duct. 
Low velocity region in free stream near injector 

* = 2.64 , H s « 8.7 MJ/kg . M ■ 3.6 . T t « 300 K 
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Pig. 18 Local equivalence ratio contours In a constant area duct 
♦ = 2.64, H s » 8.7 MJ/kg . M ■ 3.6 , T t = 300 K 




s 



Fig. 19 Axial variation of combustion and mixing efficiencies 
* = 2.64 , H = 8.7 MJ/kg , M = 3.6 , T = 300 K 
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Fig. 20 Axial variation of pressure in a constant area duct 

<t> * 2.64 . H s * 8.7 MJ/kg , U m 3.6 , T t « 1800 K 
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Pig. 22 Local equivalence ratio contours In a constant area duct 
♦ = 2.84 . H s = 8.7 MJ/kg . M = 3.6 . T^= 1800 K 
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Fig. 23 Axial variation of combustion and mixing efficiencies 
<J> * 2.64 , H s = 8.7 MJ/kg . 11 * 3.6 , T t = 1800 K 
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Fig* 24 Axial variation of pressure 0 *16* 
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Fig. 25 Axial variation of pressure e *16* 

D 


$ = 1.43, 4.2 MJ/kg . M ■ 3.5 . T t « 1800 K 


ORIGINAL PAGE IS 
Of POOR QUALITY 


89 


COMBUSTION OF HYDROGEN AND HYDROGEN/ S I LANE MIXTURES 
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SECTION 1 


COMBUSTION OF HYDROGEN/AIR MIXTURES 


This section examines hydrogen/air reactions in supersonic and hypersonic flows. Ignition 
limits have been found by experiment for a few cases (1) and these can be compared with com- 
puter simulations using a 1-D premixed finite rate chemical kinetics program (2) presented in 
this section. The number of reactions used in the simulation of hydrogen combusnon is 
reduced from 60 to 17 reactions by eliminating groups of reactions that do not appear to have a 
great effect on the combustion process at the experimental conditions. This is done in order to 

reduce computation time. 

1.1 Reactions of Hydrogen in Air 

Analysis of reactions of hydrogen in air can be initially investigated by considering only 
the oxidation of hydrogen gas into water. Using mass, momentum, and energy conservation 
equations, together with heats of reaction, final state properties can be easily calculated. These 
can be used for a quantitative comparison with experimental results of supersonic reacnng 
hydrogen/air mixtures. This approach does not take into account viscous effects, heat transfer to 
the walls, or mixing but is a useful tool in understanding the combustion process as well as giv- 
ing a reasonable approximation to final state conditions. 

A more rigorous approach, however, is needed when investigating the ignition- 
combustion limits of such a mixture. Ignition delay times become significant in very high 
speed flows where reaction speeds are critical to a scramjet’s performance. Therefore, a chemi- 
cal kinetics approach is needed to simulate the burning process. The combustion process is 
controlled by the chemical kinetics of a series of reactions whose rates can usually be found in 
the literature. For hydrogen based reactions, these rates are usually well documented. The 
difficulty lies in the choice of reactions which sufficiently model the actual combustion process. 
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A, temperatures above 1000 K, the hydrogen/air system can be closely .presented by the reac 
dons: 


H + 0 2 ^ OH + O 
0 + H 2^ 0H + H 
OH + H 2 ^ H 2 0 + H 

h 2 + m^h + h+m 

H 2 O + M ^ H + OH .+ M 

o 2 + m^o + o+m 

OH + M^O + H + M 

However, if roost of the major chemical species which comprise dr are considered, then 
them am at leas, sixty chemical reactions which can be used in the simulation of hydrogen/air 
combustion. These include the reactions that involve the gases N a and CO, in conjunction with 

O: 


By increasing the number of known 


expected to approximate more closely the combustion 


reactions in the analysis, the simulation would be 


process. However, some groups of reac- 


oons involving a single chemical specie play only a very small pan in die ignition process and 
so could be ignored to save computer time. 

In continuation of work earned out by Morgan (3) an optimal maction scheme was found 
for die combustion of hydrogen/air mixtures using a 1-D chemical kinetics analysis a, high 
temperatures. Groups of equations wem combined differently until die leas, number of .ac- 
tions were able ,0 produce results that were comparable with the sixty maction scheme. 


In all, die optimal combination comprises 17 reactions which include die seven reactions 
Stilled above together with die HO a , NO, and HNO groups of reactions. I, is noied that carbon 
species play only a very small role in hydrogen/air combustion because of the low concentre- 
tion of carbon dioxide in air. The 17 reactions together with reaction rate coefficients arc found 
in Appendix A together with the remaining reactions dtat make up the 60 reaction scheme. 
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1.2 Effect of Free Radicals on Combustion 

Previous work on the effect of ftee radicals has been done by Carson (4). Using a finite 
rate chemical kinetics program, Carson was able to demonstrate that small amounts of atomic 
oxygen, atomic hydrogen, and hydroxyl radicals reduced induction time , defined as the time to 
5% of the temperature rise from the initial mixture temperature to the calculated equilibrium 
temperature. It was also shown that hydroperoxyl (H0 2 ) plays a significant role in the early 
part of the ignition process - the inflection in its rate of formation could be used to define 
induction time. The reaction time, defined as the time between 5% and 95% of the temperature 
was not affected by the initial amount of free radicals. All of the radicals were more effective in 
reducing induction rime at low initial temperatures. For all mixtures and states, atomic oxygen 
was most effective in reducing induction time at low initial temperatures. In order to optimize 
the combustion chamber length or permit operation in a pressure or temperature limited 
regime, it was implied that additives that yield free radicals could be used. 

For work presented in this section, simulations were run at an initial nominal Mach 
number of 3.5 and static pressure of 160 kPa to coincide with experimental conditions dis- 
cussed in (1). Using a One-Dimensional Non-Equilibrium Nozzle Flow (NENZF) program (5), 
the percentage dissociation of molecular oxygen, a, was calculated for the conditions found at 
the exit of the shock tunnel nozzle. Because experimental results are presented as 
pressure/distance profiles, simulations are presented in the same way for comparison of relative 
ignition distances and relative pressure rises. Because mixture velocities stay roughly constant 
up until the point of ignition, ignition distances are directly related to ignition delay times for 
runs with identical initial conditions. 

Figure 1.1 shows a pressure-distance profile for reacting flow over 30 cm for a premixed 
stoichiometric hydrogen air mixture initially at Mach 3.5 (nominal), a temperature of 1100 K, 
and pressure of 160 kPa. These conditions correspond to an airflow stagnation enthalpy of 4.2 
MJ/kg or a flight speed of 2.9 km/s. For this case, there was no assumed initial dissociation of 
oxygen. It can be seen that, a ten reaction scheme has the longest ignition delay time. The 17, 
22, and 28 reaction schemes’ pressure profiles lay on top of each other at this condition and the 
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60 reaction scheme produces the fastest ignition. 

At these conditions, NENZF calculates an oxygen dissociation of 0.18% when the shock 
tunnel T3 is run. Figure 1.2 illustrates the effect this fraction of dissociated oxygen has on the 
simulation. The 7 and 10 reaction schemes show a marginally longer ignition delay time over 
the 17,22,28 and 60 reaction schemes. All reaction schemes show a significantly reduced 
induction time over the cases when there is no initial free radical oxygen present. 

When the temperature is reduced to 900 K (corresponding to an airflow stagnation 

enthalpy of 3.43 MJ/kg),a run over 30 cm produces no combustion when there is no initial free 

radical oxygen present (Fig. 1.3). When an initial concentration of 0.13% atomic oxygen is 

included, the same that would be present in a shock tunnel run at those conditions, combustion 

occurs within the first 10 cm as shown in Figure 1.4. Again the 7 and 10 reaction schemes fall 

upon one another and the 17,22,28 and 60 reaction schemes yield essentially the same result at 
a shorter distance. 

When the temperature is further reduced to 700K, no combustion is observed in the simu- 
lations for both cases when there is no atomic oxygen included in the analysis and with an a of 
0.06% as predicted by NENZF for a shock tunnel run at these conditions. This confirms experi- 
mental observations of no combustion at this condition, which is discussed in (1). 

The increase in ignition distance for the lower temperature case of 900K over the 1 100K 
case when a small amount of free radical oxygen is added confirms the observations by Carson 
(4). The simulations also demonstrate that the 17 reaction scheme is the lowest number of reac- 
tions that most follows the 60 reaction scheme, particularly when there is some free radical 
oxygen present. Because shock tunnels produce high speed air flows with traces of free radi- 
cals, simulations of H 2 combustion suggest that early ignition at lower temperatures in the 
shock tunnel is probably due to the presence of free radical oxygen. 

However, at higher enthalpy (hence higher temperature) conditions, the presence of larger 
amounts of free radical oxygen does not appear to have an appreciable effect on the ignition 
distance. This is illustrated in pressure-distance profiles plotted for simulations run at condi- 
tions corresponding to a stagnation enthalpy of 6.1 MJ/kg with initial pressure of 160 kPa and 
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temperature of 1700 K. Figure 1.5 shows the results for no initial atomic oxygen included in 
the analysis. The 60 reaction scheme gives marginally faster ignition over the 28,22 and 17 
reaction schemes with the 10 and 7 reaction schemes giving the largest ignition distance. At 
this enthalpy condition, NENZF calculates an oxygen dissociation of around 2.0% . When this 
is included in the combustion simulation (Fig. 1.6), the ignition is only marginally faster than 
when no atomic oxygen is included in the analysis. This suggests that if the reaction scheme is 
valid, then freestream radical oxygen production by the shock tunnel at high enthalpies does 
not have a significant effect on the ignition distance at this pressure. 

A numerical simulation of hydrogen combustion was carried out at the conditions 
corresponding with the use of the hypersonic nozzle at stagnation enthalpies of 4.2, 6.1, and 8.7 
MJ/kg and this is shown as plots of pressure against distance in Fig. 1.7. 

At 4.2 MJ/kg ( Tj =650 K ), no hydrogen combustion is observed and this is in agreement 
with experiment shown in (1). 

At 6.1 MJ/kg ( Tj = 1000 K ), the simulation predicts that ignition occurs between 30 and 
40 cm downstream from the point of injection, but hydrogen does not appear to bum at all in 
the experiment. At this condition, the ignition distance appears to be nearing the length of the 
experimental model. Experiments using a longer model would be useful to check that the igni- 
tion was merely delayed rather than completely quenched for this case. 

At 8.7 MJ/kg ( Tj = 1500 K ), the simulation predicts almost immediate combustion with 
pressure increases of between 5 and 10 kPa ( 0.05 and 0.1 atm. ). This agrees reasonably well 
with experiment shown in (1) with small pressure rises of around 5 kPa above fuel-off levels. 
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Fig. 1.3 Hydrogen Combustion Simulation. 

Effect of Varying Reactions, Mach 3.5 nominal. 

H s = 3.43 MJ/kg, Pj = 160 kPa, Tj = 900 K, <p =1.0, a = 0.00 % 
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Fig. 1.4 Hydrogen Combustion Simulation. 

Effect of Varying Reactions, Mach 3.5 nominal. 

H s = 4.2 MJ/kg, p! = 160 kPa, Ti = 900 K, 0 =1.0, a = 0 13 % 
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Fig. 1.5 Hydrogen Combustion Simulation. 

Effect of Varying Reactions, Mach 3.5 nominal. 

H s = 6.1 MJ/kg, Pj = 160 kPa, Tj = 1700 K, 0 =1.0, a = 0.00 % 



1 

1.0 i 1 j 1 1 1 I 

0 5 10 15 20 25 30 

X. cm 


60 reactions 
28 reactions 
22 reactions 
17 reactions 
10 reactions 
7 reactions 


ORIGINAL PAGE IS 
OF POOR QUALfTY 


Fig. 1.6 Hydrogen Combustion Simulation. 

Effect of Varying Reactions, Mach 3.5 nominal. 

H s = 6.1 MJ/kg. Pj = 160 kPa, Tj = 1700 K, 0 =1.0, a = 2.0 % 
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Fig. 1.7 Hydrogen Combustion Simulation. 
Effect of Varying Enthalpies, Mach 5 nominal. 
P, = 20 kPa, 0 =1-0 
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SECTION 2 


modelling the silane ignition process 


™s SeCIi ° n deSCTib ' s the r “ cdon scheroe “sed for modelling the silane/hydrogen igm- 
o°n process together with underlying assumptions. Computer simulations use a 1-D chemical 
kinetics program (2). Methods for determining the themtochemical properties of chemical 

sp^ies produced during the silme oxidation process as weU as chemical rate coefficients are 
also discussed. 

2.1 An Analogy to Methane Oxidation . The Ja c hi mowski Model. 

A chemical reaction mechanism was assembled by Beach (6) using a direct analog, with 

the methane oxidation reaction mechanism. The mechanism was tested by comparing the cal- 

culated consumption of silane with experimental results from a low temperature (500-700K) 

isothermal flow reactor, Jachimowski (7) refined this mechanism by comparing the observed 

behaviour determined in shock tube studies with tha, predicted by the mechanism. I, is tins 

work by Jachimowski on which the modelling of the combustion process in tins section is 
based. 

Basic features of the high temperature oxidation of silane were reken from die methane 
oxidation mechanism since Little chemical kinetic information is known about the silane oxida- 
tion process. The methane oxidation process has been studied for many years and the essential 
reaction paths, intermediates, and products are relatively well known. 

The silane molecule has a structure similar to that of the methane molecule, and the silane 
oxidation products ( SiO. S,0 2 ) are similar to die methane oxidation products ( CO, C0 2 ,. 
Jachimowski argued that it therefore seemed reasonable to assume that similar reactions and 
reaction intermediates would occur. It was recognised that silicon and carbon bonds were not 
identical therefore the actual reactivity and nature of the intermediate species may be different 
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from the methane scheme. However, it was the objective to develop a silane mechanism which 
contained enough detail to describe the essential features of the oxidation process rather than 
exactly model the process. 

The silane reaction mechanism used is shown in Fig. 2.1. This together with the 17 reac- 
tion hydrogen scheme outlined in section 2 as being optimum were used for any subsequent 
modelling work in this section. Rate coefficients are well documented for the current hydrogen 
schemes. However, little is known about the silane oxidation reactions. Rate coefficients 
estimated by Jachimowski using the carbon based analogy are therefore used for the work 
presented in this section. 


■Reaction 


SiH 4 + siH 2 + H 2 
SiH 4 + 0 2 -*• SiH 3 + H0 2 
SiH 4 + H0 2 -*■ SiH 3 ♦ H 2 0 2 
H + SlH 4 ■* H 2 + SiH-j 
0 + SiH 4 ♦ OH + SiHj 
OH + SiH 4 ■* H 2 0 + SiH 3 
H + SiHj -*• SiH 2 + H 2 
0 + SiH 3 -» SiH 2 0 + H 
OH + Si H j -*• SiH 2 0 + H 2 
SiH 3 + 0 2 ♦ SiH 2 0 + OH 
SiH 2 + 0 2 ♦ HSiO + OH 
H -*• SiH 2 0 ♦ H 2 + HSiO 
0 ♦ Si HjO •* OH + HSiO 
OH + SiH 2 0 + HjO + HSiO 
H + HSiO -*• H 2 + SiO 

0 + HSiO -*• OH + SiO 

OH + HS^O -» H 2 0 + SiO 
HSiO + ~M -* H + SiO + M 
HSiO + 0 2 -*■ SiO + H0 2 
SiH 2 0 + H0 2 ■* HSiO + H 2 0 2 
SiO + 0 + M -*■ Si0 2 + M 
Si O + OH -» Si0 2 + H 

Si 0 * o 2 ■» Si0 2 + 0 
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Fig. 2.1. Silane Reaction Mechanism 


2.2 Thermochemical Properties of Intermediate Species 

As indicated in (2), thermochemical properties for the program are of the form of a poly- 
nomial with seven constants describing the properties specific heat, enthalpy and entropy. The 
first five constants Aj • • • A 5 describe the specific heats while A 6 and A 7 describe standard 
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enthalpies and entropies respectively (divided by the uni 


A 6 and A 7 take the form: 


universal gas constant). In simple terms. 


A AHjJiL 


where 


-( A,T+ ^T 2 + ^ + 


>r=; 


298K 


( 1 ) 


^298 KAB) = AH 2 ° 98 ,<A) + AH 2 ° 98 KB) - D a _ 
D a _ b refers to bond energy and heats of formation 
An example of equation 3 would be 


'a-b (2 ) 

are found from assigned reference elements. 


AH 298 ,<SiH 2 0) - AH 2 ° 98 KHSiO) + AH^ g ,<H) - D H _] 
For the constant relating to standard entropy 


HSiO 


and 


c 0 

onog a * 

A 7 = ' (A > ,nT + A 2 T + ~T 2 + ~ T 3 + & 


T 4 ) 


(3) 


(4) 


S 2°98(AB) * S 2 ° 98 (A) + S& 8 (B) - AS 
where AS is the standard entropy change for the reaction 

AB -» A + B 

<n the Jachimowski model, four species SiH., SiH. SiO and Si0 2 hava documented dtar- 
mochemical ptopanias (8). The .henncchemical ptoparuas for tha othar silicon conuining 

s7ri ^ ^ "" ~ “ — - - - specie 

* S ' % S,H! °' “ d $i0 "" “ » - •*- • ** of dtair analogous carbon con- 

U 7 8 SPCC ‘“' ***** f “ SiH > - were calculated using known values 0 f 

«n alp, as of formanon and known bond anargias while stands* anthalpias for SiH 2 0 and SiO 
were calculated using known v*ues of enthalpies o, formation and ^ cncrgles h 

cerbon containing species. Standard enuopies for SiH, SiH, SiH 2 0. and SiO ware calculated 
urtng known standard enoopies and enbopy changes for analogous cart™ containing species 
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2.3 Estimation of Rate Coefficients 

The rate coefficient comprises two characteristic values for each reaction, a pre- 
exponential factor A, and an activation energy E. In the Jachimowski model, 5 reactions are 
documented and referenced with both values. The rest of the equations use the methane anal- 
ogy, drawing on either one or both of the pre-exponential factor and activation energy making 
up the rate equation. In five cases, activation energies are calculated from bond energies or 
changes in enthalpy while the rest are set equal to values found in the analogous methane reac- 
tions. 

Plots shown in Fig. 2.2 reproduced from reference (7) show the results of computer 
simulations compared with experimental shock tube results. Calculations using original 
estimated rate coefficients by Beach (6) give a much longer ignition delay time than that which 
was observed in experiment. The model was refined by Jachimowski (7) by carrying out a 
sensitivity analysis to determine which reactions had the largest effect on the ignition delay 
times. Rate coefficients were varied between 0.1 and 10 times their original value shown in Fig. 
2 . 1 . 

Three reactions appeared to have the greatest effect on the ignition delay time when their 
rate coefficients were adjusted. These were 


SiH 4 — ► SiH 2 + H 2 

(5) 

SiHj + 0 2 — > SiH 2 0 + OH 

(6) 

SiO + OH -4 Si0 2 + H 

(7) 


Each rate coefficient of each of these reactions was adjusted separately until calculated 
results agreed with experimental results. Two of the reactions yielded rate coefficients with 
unreasonably large values. Consequently it was decided that the rate coefficient should be 
adjusted for the equation which had the largest influence on calculated ignition delay times. 
This was for the reaction 

SiH 3 + 0 2 ->SiH 2 0 + OH 
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(a) 2 percent SiH^, 8 percent Hj, 4 percent 0 2 » and 
86 percent mixture; Pressure « 1.25 atm. 



(b) 1.68 percent SiH^, 6.72 percent H 2# 6.74 percent 0-, 

84.86 percent N 2 mixture; Pressure ■ 1.35 atm. 2 

Fig. 2.2. Comparison between Calculated and Experimental Ignition 


Delay Times. ' 106 
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•'»*«* T. ■ «« the ignition delay time as approximately ^ 

d ‘ S “ n “ ’° i ‘" ib0 " f " * 20% silane/hydrogen mix, urn. „ .round 70 w 

ms is much less than is predicied by - extiapolated version of Fig. 2.3 (a,. This suggest she, 

dunn, experimems. additional tew other than initial temperature and pressure of the frees 
"«m air come into play. These couid be boundary ,ayer heating on the central injection sou, 
and/or free rad, cal oxygen production in the shock tunnel which would both have an accelerate 
ing effect on the ignition process. These phenomena are discussed in section 3. 

" “ " 0ttd fr ° n ’ R8S 2 3a “ d 2 3b *»• increasing concenmttion of the silane in hydnv 
6 s a greater effect (i.e. a greater sensitization) in the stoichiometric mixtures than in the 
mixture of lower equivalence ratio. Also, ignition delay times decrease rapidly from the pure 
y g value as the amount of silane is increased. Jachitnowski noted that whereas the 
Beach mechanism (6) predicts that the same ignition delay times ate obtained for temperatures 
above 900K and for all percentages of added silane, the modified mechanism predicts tin, 
ignition delay times decrease as the silane concentration increased over the range of tempera- 
«« examined. This is supported by experimental results using the hypersonic nozzle a, 6 1 
MJ/kg (approx 1000K intake temperature, indicate in Fig. 2.4 where there is indeed a 
decrease in ignition delay rime with an increase in silane concentration. 

A calculated time histtny reproduced from (7) of silane, molecular hydrogen, and 
hydroxyl radical concentrations and mixture temperature for the stoichiometric 2% and 20% 
stlane/hydrogen mixtures a. an initial temperature of 800 K and constant pressure of 1 am, is 
shown in Fig 2.5. Jachitnowski notes that for the 2% mixture tha, the silane is consumed 
before any significant amount of hydrogen begins to react The hydrogen is then aided by the 
Urge amounts of free radicals H. O. and OH which are produced by die ignition of silane. As 
stlane concentration is increased, temperature increase due to die ignition of silane also contri- 
butes to die enhanced oxidation of die hydrogen. This diennal effect becomes mom imponant 
for sdane concentrations greater than 10%. A, 20% concentration die combined free radical and 
thermal effect is evident. After most of die silane has been consumed, and prior to significant 
oxidation of the hydrogen, the temperature has increased almost 200 K. The increase in tern- 
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P ,kP« 



■ 2.57. SiH k /H, injection, 4 * 1.08 
A 207. SiH^/H, injection, 0 = 0.90 



9 



perature and presence of free radicals both contribute to the increased rate of hydrogen oxida- 


tion. 



Infti^y ai^800 ^ istories of Sclectcd Species and Temperature During Ignition at 1 Atm and 


As well as enhancement of the ignition process through generation of free radicals, silane 
also aids ignition to some extent by scavenging the H0 2 radicals through the reaction sequence: 

SiH 4 + HO 2 — » H 2 O 2 + SiHj (g) 

H 2 O 2 + M — » OH + OH + M (9) 

At temperatures above 800 K, free radical generation and the thermal effect are the dominant 
enhancement mechanisms while below 800 K, the scavenging is probably more important (7). 

Kinetic reaction mechanisms have a hierarchical structure with mechanisms for complex 
fuels built up on sub-mechanisms for simple fuel molecules (10). The base for all hydrocarbon 
oxidation is the submechanism for the oxidation of hydrogen and carbon dioxide. The most 
important of these reactions are those that consume H 2 and CO using OH (1 1). A reasonable 
assumption that could be made for silane oxidation based on the methane analogy is that the 
important reactions would be those that consume H 2 and SiO using OH. These would be : 
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h 2 +oh->h 2 o+h 


( 10 ) 


SiO + OH — ► Si0 2 + H (11) 

together with the dominant chain branching reaction 

H + 0 2 -» O + OH (12) 

and its primary competitor for H atoms 

H + 0 2 + M->H0 2 + M (13) 

Reactions competing with H 2 and SiO for OH should inhibit or retard oxidation of these 
species as well as the heat release associated with the production of water and silicon dioxide, 
and reactions competing with reaction 12 for H atoms (reaction 13) should reduce the rate of 
chain branching and overall rate of combustion. Alternatively, reactions which produce addi- 
tional H atoms which can then react with molecular oxygen by reaction 12 should accelerate 
the overall rate of combustion. 

It will be seen, however, in section 3.5 that according to a numerical chemical kinetics 
analysis, reaction equations 10 and 11 do not play any significant part in the ignition process 
for particular initial conditions. In fact, hydroxyl radicals are important in producing intermedi- 
ate silicon containing species but not silicon dioxide. Heat release is achieved by a series of 
reactions culminating in the direct attack of silicon monoxide on molecular oxygen and then by 
the ignition of hydrogen. At lower temperatures, hydrogen stays essentially "inert" and does 
take pan to any significant degree in the ignition process. Direct modelling on a hydrocarbon 
analogy, therefore, has its drawbacks. Conclusions which can be drawn from modelling work 
concerning methane, for example, cannot be always applied directly to silane combustion even 
though the silane reaction scheme is based on an equivalent methane reaction scheme. 

The oxidation of silane takes place through the gradual dismantling of the fuel molecule 
into progressively smaller fragments. The first major step is the abstraction of H atoms from 
the parent fuel molecule. Once a viable radical pool has been established, most of the H atom 
abstraction is accomplished by reactions between small radical and atomic species and the fuel 
molecules. The most important radical species are H, O, OH, H0 2 , and SiH 3 . From previous 
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hydrocarbon modelling work, H atom abstractions by H and OH are dominant for fuel rich and 
near stoichiometric mixtures while reactions with OH and to a much lesser extent, O atoms are 
most important for fuel lean mixtures (11). Work in section 3.5 tends to support the trend of H 
and especially OH radicals being responsible for most hydrogen abstraction from the parent 
silane molecule when there is a stoichiometric fuel/air mixture. In this particular manner, silane- 
tends to behave like its hydrocarbon counterpart 

2.5 An Alternative Method for Approximating Thermochemical Properties 

The Jachimowski model approximates for 4 of the silicon species involved in the silane 
oxidation process ( SiH 3 , SiH 2 , SiH 2 0, HSiO ) that the heat capacities are the same as those for 
the analogous carbon containing species. An indication of the possible difference between 
assumed values and actual values can be demonstrated by comparing specific heats of say SiH 4 
and CH 4 . For example, at 1000 K, the specific heat of SiH 4 is 20.2 cal/mole K while that of 
CH 4 is 17.4 cal/mole K which represents an eiror of 14% if the analogous carbon containing 
specie was used in this case. It was thought that an independent method of determining heat 
capacities could be investigated to check the effects on simulations. 

Using the method of Bennewitz, Rossner and Dobratz (12), molecular specific heats can 
be determined approximately by considering contributions due to translational and rotational 
energies together with vibrational energies expressed as functions of temperature. Using the 
principle of equipartition, translational and rotational contributions to C p are 3R. Vibrational 
contributions depend on bond frequencies in stretching or bending. Average stretching or 
bending frequencies are expressed as wave, numbers (frequency divided by the speed of light), 
where to, is the wave number for stretching and co 5 is the wave number for bending. Vibra- 
tional contributions to heat capacities are related to these wave numbers independent of the 

bond type and take the form of functions dependent on temperature. These contributions are 
then multiplied for multiple bonds. 

The average stretching vibration wave number for the C-H bond is o\, = 2920 cm -1 (12). 
For Si-H. co, = 2190 cm -1 (8). This represents a difference of 730cm-' in wave number 
between analogous species. For the sulphur bond S-H, tn, = 2570 cm-' (12) which represents a 
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smaller difference of 380 cm -1 . It would be expected therefore that using S-H bonds instead of 
C-H bonds would give a better approximation to the specific heat of the silicon containing 
specie. This is indeed the case and at 1000 K, Cp calculated for SiH 4 using this "sulphur anal- 
ogy" is 19.5 cal/mole K compared with the documented value of 20.2 cal/mole K, an error of 
3%. 

The model is even more accurate for the documented specie Si=0. For Si=0 bonds, 
toy = 1295 cm -1 (8), while for S=0 bonds o\,= 1250 cm -1 (12), representing a difference in 
wave number of only 45 cm -1 . Calculated specific heat using the "sulphur analogy" for SiO at 
1000 K is 8.56 cal/ mole K compared with the documented value of 8.54 cal/mole K, an error 
of only 0.2%. This compares with an error of 7% when using the analogous carbon containing 
specie CO. 

It was on this basis that specific heats were calculated for the undocumented silicon 
species using a sulphur bond analogy. Standard heats of formation were calculated using the 
method of Franklin, Verma and Doraiswarmy (12), and absolute entropies using the same 
method as Jachimowski except that entropy changes for reactions were approximated by sums 
of known entropy changes involved in the breaking of single bonds rather than using the carbon 
reaction analogy. 

The modified thermochemical properties were used for numerical simulations in conjunc- 
tion with reaction rate coefficients suggested by Jachimowski. The modified properties gave 
similar results for the Jachimowski model for several cases, with slightly longer ignition delay 
times but similar final pressures. Computation time was increased gready and reasons for this 
are not known. Overall, it could be concluded that small errors in calculating thermochemical 
properties have a far less effect on simulation results than choice of reaction rate constants. 
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SECTION 3 


SIMULATION OF THE SILANE COMBUSTION PROCESS 


Usmg the Jachimowski model described in section 2 and the chemical kinetics program 
(2), computer simulations have been carried out for conditions corresponding to experiments 
using the hypersonic nozzle. Experimental fuel off duct pressures varied between 13-22 kPa 
over the enthalpy range considered but simulations were consistently run with an initial pres- 
sure of 20 kPa. Effects of freestream free radical oxygen concentration, concentration of silane 
in hydrogen, and equivalence ratio are discussed and results compared with experiment. A sen- 
sitivity analysis of reaction rates is carried out at a condition where results of experiment 
correspond well with experimental results. An examination of net conversion rates at this same 
condition indicates important reactions at various stages of the ignition process. 

3.1 Effect of Freestream Free Radical Oxygen Concentration 

All work presented in this section concerns combustion of a 20% silane^ydrogen mixture 
at an equivalence ratio of one. Fig. 3.1 shows that at 1500 K or a freestream stagnation 
enthalpy of 8.7 MJ/kg, increasing the initial concentration of freestream free radical oxygen has 
little or no effect on the ignition distances. Results from this simulation compare favourably 
with the experimental results for a 20% mixture as shown in Fig. 3.2. Note that for computer 
simulations the point X=0 refers to the point of injection, rather than the position of the leading 
edge. The experimental points appear scattered, but at these high initial temperatures, pressure 
rises due to combustion are relatively low, so any disturbances within the model would look 
significant against the the combustion profile. It can be seen, however, that there is indeed rapid 
combustion that comes up to roughly the level predicted by the computer simulation. It is sug- 
gested therefore that at this high enthalpy, the high dissociation expected in the shock tunnel 
(a = 10%) has no noticeable effect on the ignition delay time. 
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Fig. 3.1 20% Silane/Hydrogen Combustion Simulation. 

Effect of Concentration of Free Radical Oxygen. Mach 5 no min al 
H s = 8.7 MJ/kg, Pi = 20 kPa, Tj = 1500 K, cf> =1.0 
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Fig. 3.2 20% Silane/Hydrogen Combustion. 

Comparison of Simulation with Experiment. Mac 
H s = 8.7 MJ/kg, pj = 20 kPa, T, = 1500 K, $ =1.( 


5 nominal, 
a = 10% 
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At an initial temperature of 1000 Kora freestream stagnation enthalpy of 6.1 MJ/kg. 
simulations show the ignition distance decreases from just over 20 cm to 5 cm as the initial a is 
increased from 0.0% to 0.8% (Fig. 3.3). At this condition, the shock tunnel produces flows with 
a = 1.5%. According to the simulation, this would bring the rapid pressure rise to around 20 
cm closer to the point of injection than when compared to a flow with no oxygen dissociation. 
Fig. 3.4 shows that in the experiment, results agree favourably with the simulation run with 
with an initial free radical oxygen concentration of 1.5%. It is interesting to note that in experi- 
ment, the pressure falls off with distance down the duct. 

Fig. 3.5 shows that at an initial temperature of 650 K ( H s = 4.2 MJ/kg ), even small 
amounts of free radical oxygen have a very significant effect on the ignition distances predicted 
by the computer simulation. At this enthalpy a = 0.18% is expected in the shock tunnel. It can 
be seen by comparing the results for a 20% silane mixture shown in Fig. 3.6 with the simula- 
tion for a = 0.18% that ignition distance is about the same. At this temperature as at the higher 
temperature of 1000 K, a pressure drop is observed as the gas progresses down the duct This 
trend appears to be repeated for most cases with silane injection, particularly at the higher con- 
centrations. A possible cause of this may be that after ignition and burning of the 
silane/hydrogen mixture is complete, heat transfer to the walls of the model may cause conden- 
sation of Si0 2 . This conversion from gas to solid would tend to lower pressure with distance 
more drastically than the effect due to a drop in temperature alone, as found in say, pure hydro- 
gen combustion. Results for a 20% silane/hydrogen mixture at nominal Mach numbers of 
Mach 3.5 and 4.5 show a similar trend as shown in (1,9). It should be noted that for hydrogen 
combustion, where water remains as a vapour at relatively low temperatures, that there is little 
pressure fall off as found in (1). 

At a lower initial temperature of 540 K ( H s = 3.43 MJ/kg ), ignition distances are dramat- 
ically increased for the lower radical concentrations, as predicted by the computer simulation 
(Fig. 3.7). However, Fig. 3.8 shows that in experiments with o = 0.13%, the 20% silane mix- 
ture ignites about 15 cm from the point of injection. This compares with about 65 cm for the 
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Fig. 3.3 20% Silane/Hydrogen Combustion Simulation. 

Effect of Concentration of Free Radical Oxygen. Mach 5 nominal. 
H s = 6.1 MJ/kg, Pi = 20 kPa, Tj = 1000 K, <t> =1 .0 
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Fig. 3.4 20% Silane/Hydrogen Combustion. 

Comparison of Simulation with Experiment. Mach 5 nominal. 
H s = 6.1 MJ/kg, Pj = 20 kPa, T, = 1000 K, $ =1.0, a = 1.5% 
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Fig. 3.5 20% Silane/Hydrogen Combustion Simulation. 

Effect of Concentration of Free Radical Oxygen. Mach 5 nominal. 
H s = 4.2 MJ/kg, Pj = 20 kPa, Tj = 650 K, <J> =1 .0 
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Fig. 3.6 20% Silane/Hydrogen Combustion. 

Comparison of Simulation with E xperiment. Mach 5 nominal. 
H s - 4.2 MJ/kg, Pj = 20 kPa, T, = 650 K, $ =1 .0, a = 0. 1 8% 
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Fig. 3.7 20% Silane/Hydrogen Combustion Simulation. 

Effect of Concentration of Free Radical Oxygen. Mach 5 nominal. 
H s = 3.43 MJ/kg, P, = 20 kPa, Tj = 540 K, 0 =1 .0 
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Fig. 3.8 20% Silane/Hydrogen Combustion. 

Comparison of Simulation with Experiment. Mach 5 nominal. 
H s = 3.43 MJ/kg, P, = 20 kPa, Tj = 540 K, 0=1.0, a = 10% 



premixed gases as prediced by ft. 1-D prepare, The pressure rise ^ cxpenme „, 

is. however, only of the order of 20 kPa or 0.2 atm compared with 100 IcPa or 1.0 atm 
predicted for -full- burning of the silane/hydrogen mixture as shown i„ Rg. 3.7. This suggests 
d»t there is penial burning of the fuel je, and that ignition is caused prematurely by some other 
effect because the freestream temperature is too low. The ignition at 540 K may be due to high 
temperatures produced in the boundary layer on the injection smut. Assuming a laminar 
compressible boundary layer, i, was calculated that the maximum temperature found in the 
boundary layer would be .060 K a, this condition. This temperature together with a free radical 
oxygen concentration of 0.13 % found in the freestream would be consistent with the ignition 
distance of around 10 cm predicted in Rg. 3.2. The partial burning of the fuel jet could be 
caused by the quenching of the complete combustion process by the cooler freestream. 

It should be noted that i, was a. this inlet temperature, pressure, and Mach number that 
higher steady pressures were achieved experimentally with lower silane concentrations. I, 
appears therefore, burning was incomplete but involved both components of the fuel mixture. 
This effect is not understood and is not modelled by the chemical kinetics program in this case. 

A simulation was run (including o = 0.06% calculated for conditions found at the exit of 
the nozzle) for a lower intake temperature of 410 K ( H s -2.65 MJ/kg ) and indicated no 
combustion over 120 cm. Experiments, however, showed ignition at about 15 cm from the 
injection point (9) but final pressure levels did not reflect complete combustion. I, was 
hypothesised that preferential burning of the silane component of the fuel mixture occurs at 
these low temperatures and this is supported by experiments a. Mach 4.5 and Mach 5. The 
chemical kinetics program, however, does not appear to model this behaviour. 

Experiments suggest therefore that if the reaction scheme is valid, then freestream radical 

production at high enthalpies has no appreciable effect on the ignition of 20% silane/hydrogen 
mixtures. 


At the intermediate freestream temperature of 650 K, oxygen radical concentrations play 

“ i “ POrUm pan in the com hustion process therefore shock tunnel simulations may no, be 
representative of real flight situations. 
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At the lower intake temperatures of 540 K and 410 K, experiments show short ignition 
delay times despite the very low free radical oxygen concentrations produced by the shock tun- 
nel. This suggests that combustion is occurring because of another effect. It may be the high 
temperature produced by the boundary layer on the injection stmt. Pressure rises due to burning 
however, are lower than those expected for complete combustion of the silane/hydrogen mix- 
ture as predicted by the computer program. This suggested that there is partial burning due to 
the low freestream temperature or preferential burning of the silane component of the fuel mix- 
ture. 


3.2 Effect of Silane Concentration 

Conditions which had previously shown good agreement with experiment were chosen for 
investigating the effect of varying the concentration of silane in hydrogen. It has been shown 
that at a freestream stagnation enthalpy of 4.2 MJ/kg, the simulation compares favourably with 
experiment for a 20% mixture when free radical oxygen produced in the shock tunnel is 
included in the analysis. 

Fig. 3.9 shows that at this condition, pressure profiles indicate a change in ignition trends 
as concentration is increased. At 1%, there is minimal ignition over 80 cm. As the concentra- 
tion is increased from 2.5% to 20%, ignition distances are decreased and the ignition process is 
indicated by more rapid pressure rises. This compares favourably with experiment as shown in 
(9) where similar trends are observed, although ignition distances tend to be smaller 

3.3 Effect of Equivalence Ratio 

Equivalence ratio was varied for the same condition. Fig. 3.10 shows that ignition delay 
times do not vary much up to an equivalence ratio of about 2.0. However, at an equivalence 
ratio of 5.0, ignition delay is increased significantly. This could be due to the heat capacity of 
the excess fuel robbing the ignition process of the necessary heat. Final pressure levels are 
about the same for ratios between 1 .0 and 2.0. Below equivalence ratios of 1 .0, final pressure 
levels are sensitive to small changes in 4>. This trend was observed in expe rime nts nominally at 
Mach 4.5, where equivalence ratios were nominally 0.6. 
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Fig. 3.9 Silane/Hydrogen Combustion Simulation. 

Effect of Concentration of Silane. Mach 5 nominal. 

H s = 4.2 MJ/kg, Pj = 20 kPa, Tj = 650 K, <J> =1 .0, a = 0.18% 
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Fig. 3.10 20% Silane/Hydrogen Combustion Simulation; 

Effect of Equivalence Ratio. Mach 5 nominal 

H s = 4.2 MJ/kg, Pj = 20 kPa, T, = 650 K, 4 =1 .0, a = 0. 1 8% 
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3.4 A Sensitivity Analysis of Reaction Rates 

The ignition mechanism tends to be made up of four processes. These are: 

1. Breakdown of the silane to produce free radicals. 

2. Production of intermediate species. 

3. Ignition of silane to produce heat. 

4. Subsequent ignition of remaining hydrogen. 

A sensitivity analysis of reaction rates for reactions involving silicon containing species is 
presented in this section. To reduce complexity, results are presented with reference only to the 
23 reactions concerning silane oxidation. This is a reasonable approach because with the rela- 
tively high concentration of 20% silane in hydrogen considered for this case, the reactions 
involving silane are more or less complete before the hydrogen oxidation begins (section 2.4). 
The reactions involving silicon therefore tend to be inter-related with one another to a greater 
extent than with the 17 reactions which describe the hydrogen oxidation. 

Rate coefficients of the silane reaction scheme shown in Fig. 2.1 were successively varied 
by factors of 2 and 0.5 for the simulation that gave best agreement with the experimental 
results. This was at the condition discussed in section 3.2 with H s = 4.2 MJ/kg, Tj = 650 K and 
Pi = 20kPa with the free radical oxygen included in the analysis. The sensitivity analysis was 
based on atomic oxygen mass fractions computed for the experimental conditions at distances 
along the duct. Atomic oxygen was chosen as the parameter for observation because of its 
apparent importance in the ignition process at this condition. Plots of mass fraction against 
length corresponding to distance from injection are shown in Figs. 3.11 (a) to (i), with rate 
coefficients kl to k23 corresponding to the reactions R1 to R23 shown in the accompanying 
table. The result for this condition using original rate coefficients is shown by a continuous line. 
The importance of each reaction is signified by the deviation from the original mass fraction 
profile when the rate coefficient corresponding to that reaction is halved or doubled. Points 
corresponding to the most significant reactions have been connected with broken lines. 

From pressure-distance profiles which are not presented here, it was found that ignition 
(as signified by a rapid rise in pressure) consistently took place when the pool of free radical 
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Fig. 3.1 1 20% Silane/Hydrogen Combustion Simulation 

H ^ Ra,cs - Mach 5 nominal. 

H s 4.2 MJ/kg, P, _ 20 kPa, T, = 650 K, <> =1 .0, a = 0. 1 8% 
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Fig. 3.1 1 20% Silane/I lydrogen Combustion Simujation. 

Sensitivity Analysis of Reaction Rates. Mach 5 nominal. 

H s = 4.2 MJ/kg. P, = 20 kPa, T, = 650 K, <> =1 .0, a = 0. 1 8% 



Rl 

R2 

R3 

R4 

R5 

R6 

R7 

R8 

R9 

RIO 

Rll 

R12 

R13 

R14 

R15 

R16 

R17 

R18 

R19 

R20 
R21 
R22 
R2 3 


Reaction 

SiH 4 -► SiH 2 + H 2 
siH 4 + 0 2 -*■ SiH 3 + H0 2 
SiH 4 + H0 2 -► SiH 3 +' H-,0, 
H + SiH 4 ^ h 2 + SiH 3 
O + SiH 4 OH + SiH 3 
OH + SiH 4 -► h 2 0 + SiH 3 
H + SiH 3 ■+• SiH 2 + h 2 
0 + SiH 3 SiH 2 0 + H 
OH + SiH 3 + SiH 2 0 + H 2 
SiH 3 + 0 2 -*• SiH 2 0 + OH 
SiH 2 + 0 2 -► HSiO + OH 
H + SiH^O -► H 2 + HSiO 
0 + SiH 2 0 -► OH + HSiO 
OH + SiH 2 o -► H 2 0 + HSiO 
H + HSiO ♦ H 2 + SiO 
0 + HSiO -► OH + SiO 
OH + HSijO -► H 2 0 + SiO 
HSiO + -*• H + SiO + M 

HSiO + 0 2 -► SiO + HO 2 
SiH 2 0 + H0 2 ■+ HSiO + H 2 0 2 
SiO + 0 + M -*• Si0 2 + m * 
SiO + OH ■* Si0 2 + H 
SiO + 0 2 ■+ Si0 2 + 0 


Rate coefficient 

6.0 x 10 13 exp ( -54960/RT) 

2.0 x io 11 exp( -44000/RT) 

3.0 x i o 1 2 exp( -5600/RT) 
1.5 x io 13 exp( -2500/RT) 

4.2 x io 12 exp(-1 600/RT) 

8.4 x io 12 exp(-100/RT) 

1.5 x io 13 exp( -2500/RT) 

1.3 x io 14 exp( -2000/RT) 

5.0 x to 1 2 

6.6 x io 14 exp(-1 1400/RT) 

1.0 x io 14 exp( -3700/RT) 

3.3 x io 14 exp( -1 0500/RT) 
1.8 x io 13 exp ( -3080/RT) 

7.5 x io 12 exp( -1 70/RT) 

2.0 x io 14 

1 .0 x io 14 

1.0 x io 14 

5.0 x io 14 exp(-29000/RT) 

3.0 x io 1 2 

1.0 x 10 12 exp( -8000/RT) 

2.5 x io 15 exp ( -4370/RT) 

4.0 x to 12 exp{ -5700/RT) 

1.0 x i o 1 3 exp( -6500/RT) 


(i) 
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oxygen reached a "critical" mass fraction of roughly 3x10"* . Ignition distances can be approxi- 
mated by considering the point where the atomic oxygen mass fraction reaches this value. This 
means generally that if changing a reaction rate increased radical oxygen mass fractions above 
that of the original scheme to the "critical" value then the ignition distance was shorter than for 
the original scheme. Alternatively, if the original scheme reached this value before the 
modified scheme, then the ignition distance for the modified reaction was shorter than for the 
original scheme. 

From the results, it can be seen that there are six reactions which have the largest 
influence when rate coefficients are varied. In order of importance from highest to lowest, 

these are 

RIO...S1H3 + Oj — * SiH 2 0 + OH 
R5...0 + SiH4 -* OH + SiH 3 
R8...0 + SiH 3 -» SiH 2 0 + H 
R14...0H + SiH 2 0 — * H 2 0 + HSiO 
R23...SiO + 0 2 -» Si0 2 + O 
R19...HSiO + 0 2 -» SiO + H0 2 

It can be seen from Figs. 3.1 1(b) and 3.1 1(f) that varying the reaction rate of RIO has the 
greatest effect on the mass fraction profile. This is in agreement with the observation of Jachi- 
mowski of greatest sensitivity for this reaction as discussed in section 2.3. It should be 
remembered that it was this reaction’s rate coefficient that Jachimowski adjusted so that numer- 
ical simulations corresponded with experimental shock tube results. The other two reactions 
which Jachimowski found to be most sensitive to changes in reaction rate were: 

Rl...SiH 4 -4SiH 2 + H 2 
R22...SiO + OH — » Si0 2 + H 

These reactions were not observed to be important in this sensitivity analysis. This is 
understandable because a small concentration of free radical oxygen was included in this par- 
ticular analysis. Reactions involving atomic oxygen would therefore tend to become dominant 
1 reactions, as is indeed observed in this case with R5 and R8 being most sensitive after RIO. 



" Sh °“ M ^ " 0Kd R8 * <•— ^ when i,s rate coefficient was do „. 

bled end *«I«a,ed the process when its rate coefficient was halved. This is though, to be 

because R8 cotnpe.es with other reactions for the radicals SiH, a* o to produce the intennedi. 

«e S. Hl 0 and atomic hydrogen. During the ignition process, SiH, and O appear to be very 

unponan, in attacking oxygen and silane respectively in sections RIO and R5. and so con- 

sumpnon of these radicals through other reactions would tend to reduce the rate of combusdon. 

TTte remaining reacdons R14, R23 and R,9 appear be imponan, macrons bu, to a 
lesser degree than RIO, R5 and R8. Except for the reacdon : 

R10...SiH 3 + 0 2 — » SiH 2 0 + OH 

dtose that were observed to be sensidve for this case diffeted from Utose dta, Jachimowski 
observed, as described in secnon 2.3. It would appear therefore that choice of inidal conditions 

h “ a tegC b ' an " g °" *• resul “ a analysis. In this cam, the most imponan. 

parameter tha, affected the outcome of the analysis was the inclusion of an inidal amount of 
free radical oxygen. 

3.5 A Comparison of Net Reaction Conversion Rates 

In this section, net forward reacdons nt.es (X,) are ploned against distance (in Figs. 3.12 

W .0 3.12 (d» for die same example as discussed in sections 3.2, 3.3, and 3.4. The magnitude 

of X; gives approximately the imponance of any single reaction amongst all those occuning. I, 
is defined for different types of reactions in (2). 

From these plots, it can be seen which reactions me relatively imponan, a„ d a , which 

““ m die ignition process they occur. The rapid pressure rim signifying compete combus- 

non occurs between 25 and 27 cm downstream from die starting point for the analysis (Fig. 

3.13). Given this, i, can be men dia, them are four distinct processes which make up die igni- 
tion process. 

For die firs, sage of die ignition process, atomic oxygen is rapidly consumed by die reac 
tion 

R5...0 + SiHj -» 0H + SiH-> 
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Net Forward Reaction Rate no)e/cm3/g2/sec Net Forward Reaction Rate Mole/cm3/g2/sec 



O Reset ion i 
O Reset ion ? 
o Resctlon 3 
A Resction 4 
v Resction 5 
x Reaction 6 



o Reaction 7 
O Reaction e 
o Reaction 9 
A Reaction 10 
^ Reaction 11 
x Reaction 12 


Fig. 3.12 20% Silane/Hydrogen Combustion Simulation. 

Comparison of Net Forward Reaction Rates. Mach 5 nominal. 
H s = 4.2 MJ/kg, Pj = 20 kPa, T, = 650 K, $ =1.0, a = 0.18% 
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Net Forward Reaction Rate mo le/c» 3 /g 2 /sec Net Forward Reaction Rate Niole/ci* 3 /g 2 /tec 



O Reaction 13 
O Reaction 14 
o Reaction 15 

* Reaction 16 
v Reaction 17 

* Reaction ia 


(d) 



o Reaction 19 
© Reaction 20 
o Reaction 2i 
A Reaction 22 
v Reaction 23 


Fig. 3. 12 20% Silane/Hydrogen Combustion Simulation. 

Companson of Net Forward Reaction Rates. Mach 5 nominal. 
H s - 4.2 MJ/kg, Pj = 20 kPa, = 650 K, <J> = 1 . 0 , a = 0. 1 8% 
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and the conversion rate for this reaction drops to about zero at about 6cm. This corresponds 
with the point of minimum free radical concentration observed for this case as illustrated by the 
atomic oxygen mass fraction profiles in section 3.4 for the unadjusted reaction scheme. After 
this initial attack by the free radical oxygen on the silane, a pool of hydroxyl andSiH 3 radicals 
is established to feed the two reactions 

R6...0H + SiH4 -> H 2 0 + SiH 3 
R10...SiH 3 + 0 2 -► SiH 2 0 + OH 

This is the second stage of the ignition process where it can be seen that one reaction sus- 
tains the other by producing the radicals the other needs. R6 provides RIO with the SiH 3 to 
attack the oxygen while RIO provides R6 with the OH to attack the silane. These two reactions 
become the dominant reactions for the early part of the ignition process up until about 15 cm. 

Another important reaction ( but to a lesser extent ) during these early stages is 

R4...H + SiH 4 H 2 + SiH 3 

The atomic hydrogen which attacks the silane is produced from the reaction 

O + H 2 -> OH + H 

This is an example of the interelation of the silane oxidation process with a reaction involving 
the hydrogen component of the fuel mixture. The overall process at this stage, however, is 
dominated by the reactions R6 and RIO which involve the silane component of the fuel mix- 
ture. 

Hydrogen abstractions from the parent fuel molecule by OH and H radicals tend to be the 
dominant reactions throughout the ignition process. Abstraction by O radicals produced by the 
shock tunnel accelerates the process significantly as observed in the first stage of the ignition 
process, but veiy little atomic oxygen takes pan in hydrogen abstraction reaction R5 thereafter. 
This trend supports previous experience in hydrocarbon modelling work with stoichiometric 
fuel mixtures, as discussed in section 2.4. 

During the second stage of ignition process, free radicals O, OH, H and SiH 3 are gen- 
erated together with a quantity of the intermediate specie SiH 2 0 through the reaction RIO. 


After 15 cm, the third stage of the reaction process becomes apparent. Three interested reac- 
tions stan to become important and net conversion rates for these three reactions rise at about 
the same rate. The three reactions are 

R14...0H + SiH 2 0 -> OH + HSiO 
R19...HSiO + Oj SiO + HO 2 
R23...SiO + 0 2 -» Si0 2 + O 

It can be seen that these three reactions complete a reaction sequence that converts silane 
into silicon dioxide. The net release of heat from this process would then be responsible for the 
increased activity of almost all of the other reactions at about 20 cm. Net conversion rates then 
nse to a maximum for most reactions at about 25 cm. This distance corresponds with the 
"hydrogen ignition" point and the fourth stage of the process where the temperature is sufficient 

to cause the hydrogen to bum. This is evidenced by the rapid pressure rise between 25 and 27 
cm as illustrated in Fig. 3.6. 
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Nomenclature 

Stagnation Enthalpy, MJ/kg 

Intake Mach Number 

Intake Pressure, kPa 

Average Fuel-off Duct Pressure, kPa 

Stagnation Pressure, kPa 

Local Static Pressure, kPa 

Equivalence Ratio 

Intake Static Temperature, K 

Duct Wetted Length, mm or cm 

Net Forward Reaction Rate, mole/cmVg^s 

Percentage Dissociation of Molecular Oxygen 
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abstract. 

The results are presented of a series of 
tests on a sodel scrasjet perforsed in 

the free piston shock tunnel * T3* in the 
Australian National University. The 
experiments were planned to five 

performance data over a range of pressure 
levels corresponding to flight at 

different altitudes. For each pressure 
level tests were carried out for 
enthalpies ranging from the ignition 
limits at the low temperature end. up to 
temperatures where the dissociation of 
combustion products severely limited heat 
release. The minimua temperature at which 
combustion was possible was found to be 
highly pressure sensitive for the 

experimental conditions used. At lower 
pressures higher temperatures were needed 
to initiate combustion within the transit 
time of the combustion chamber. Close to. 
but above, the ignition limit the amount 
; o* heat release increased markedly with 
■ Pressure and also with combustion 
i chamber length. At higher enthalpies 
! combustion was not so sensitive to 
pressure rise. A finite element computer 
code was used to model the mixing and 
combustion processes. Away from ignition 
delay effects the program gives 
reasonable agreement with the measured 
heat release. At conditions where 
combustion is marginal, near the ignition 
limits, good agreement with experiment 
can only be obtained by adjusting the 
freestream oxygen dissociation fraction 
to unreasonable levels. 


For combustion of fuel in a scrasjet 
three processes are necessary, namely 
mixing of the fuel and air, heating of 
the fuel above its ignition temperature, 
and the allowance of sufficient time for 
reactions to take place. All of these 
processes are to a certain extent 

pressure dependent. An accelerating 
will experience a wide range of 
operating conditions. These include 

pressure. Mach number , and static 

temperature variations at the combustion 

♦ Seuor Uctirtr, Mbamcal ir*irmrim, Kmtor AIM* 
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chamber intake. An understanding of the 
behaviour of scrasjets over the full 
range of flight conditions is therefore 
important . 

Previous shock tunnel studies of 
mcramjets have concentrated on basic 
thrust producing and heat transfer 
mechanisms, Refs l and 2. This paper 
reports the results of experiments to 
study pressure effects in a shock tunnel, 
and presents comparisons with a currently 
available chemical kinetics and mixing 
program* Ref 3. 

The experiments were devised to give 
■Pacific impulse against temperature data 
for a range of intake static pressures. 
From this data it was possible to 

identify three broad combustion 

categories in a pressure / temperature 
plane. Fig 9. In region 1* combustion is 
not possible within the confines of the 
chosen model. In region 2 combustion is 
possible. but the heat release is 

strongly pressure and temperature 

dependent. in region 3, combustion is 
only slightly sensitive to pressure and 
^•*P* r *ture, over the pressure range 
considered . 

Somm tests were perforsed with ah 
extension block fitted to the combustion 
chamber, which allowed more time for 
mixing and combustion before expanding 
the flow in the thrust nozzle. This was 
be effective in the combustion 
category 2, above* but had no measurable 
effect at the higher temperatures 
associated with category 3. it was also 
seen to shift the border of regions i and 
2 to lower temperatures* indicating that 
the extra mixing length permitted 
combustion under conditions with a longer 
ignition delay time. 


IXPKfllMPfTAL APPARATUS 


A schematic diagram 
apparatus is shown 
piston shock tunnel . 
Australian National 
to create the 
corresponding to the 
combustion chamber. 


of the experimental 
in Fig l. The free 
T3 Ref 2, at the 
University was used 
flow conditions 
intake of a scrasjet 


j The sodel consisted of a two dimensional 
duct of intake section 25 mm by 50 sa . 

; Hydrogen fuel was injected on the 

centreline from an injector strut which 
spanned the whole width of the duct The 
leading edge of the injector was extended 
clear of the intake so that disturbances 
from the leading edge would not enter the 
; duct. 
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After injection a constant area section 
of adjustable length served as a 
combustion chamber. The flow was expanded 
by a single 15 degree Prandtl Meyer 
expansion fan propagating from the lower 
wall . and thrust was developed on a 
straight thrust surface. The intake, 
combustion chamber and thrust surface 
were instrumented with pietzotronic 
pressure transducers. Thrust was 
calculated from the static pressure 
measurements, and specific impulse was 
determined * from the difference between 
fuel on and fuel off thrust levels. 

A steady flow of hydrogen was established 
prior to starting the shock tunnel by 
means of a pulsed valve. Ref 3. which was 
triggered off the recoil of the tunnel. A 
static pressure transducer in the 
injector plenum chamber was used to give 
a measure of the hydrogen sass flow 
rate. 

A contoured nozzle of nominal Mach number 
3.5 was used to create the required test 
conditions. The model intake static 
pressure was controlled by selection of a \ 
suitable shock tube driver rupture j 
pressure. Further control on the intake j 
pressure and Mach number was given by the j 
use of adjustable wedges on the model 1 
inlet. These wedges created oblique 
shocks which reflected from the 
centreline and did not enter the model, 
thereby maintaining uniform intake flow. 
Vt should be noted that the use of the 
wedges was purely to give flexibility to 
the conditions achievable with a single 
contoured nozzle. They do not represent 
an attempt to model a scram jet intake, 
because, due to Reynolds number 
limitations, it was only possible to 
* model the combustion chamber and a short 
thrust nozzle. 


CALCULATION OF TEST COttPITK KSr 

•The flow conditions at the nozzle exit 
are estimated assuming equilibrium in 
the shock tube stagnation region behind 
the reflected shock, followed by a 
non-equilibrium expansion through the 
contoured nozzle. The non equilibrium 
expansion is computed using the computer 
code of Ref 6 . 

Frozen flow is assumed through the 
oblique shocks formed by the compression 
wedges , when used . 

A total of 44 different shock tunnel 
operating conditions were used, all of 
which could be configured with 3 
different wedge compression angles. It 


was impractical to perform a complete 
analysis of all the intake data . so 
interpolation was used for some cases. In 
Fig 2 a range of computed nozzle exit 
velocities against stagnation enthalpy is 
plotted, and in Fig 3 the curve fit used 
for the full range of conditions is 
shown. Fig 4 shows the dependence of Mach 
number with stagnation enthalpy . The Mach 
number is reduced at higher enthalpies 
due to the a reduction in the specific 
heat ratio. 

The test conditions are tabulated in 
table 1 . Equivalence ratios were 
nominally set to 1. with minor 
fluctuations between runs . 


The shock tunnel was operated in the 
undertailored mode, which was found to 
maximise the uncontaminated test time, 
which is especially important in 
combustion applications. This condition 
produces a small drop in the stagnation 
pressure during the run. In order to 
confirm that the expansion through the 
■ nozzle produced a steady Mach number , the 
static pressures in the model were 
normalised with respect to the shock 
tunnel stagnation pressure. Readings were 
taken onlv when all the normalised 
static pressures were steady. All 
computations performed on the data 
assumed steady flow conditions. 

Specific impulse was calculated by taking 
the difference between fuel on and fuel 
off thrust levels. normalised by 
stagnation pressure to account for minor 
fluctuations between nominally identical 
runs . 

In Fig 5 the effects of pressure at a 
Mach number of 3.25 are shown for a short 
combustion chamber. A 5 degree angle was 
used on the compression wedges to create 
this condition. 

A striking feature of this figure is the 
dependence of ignition temperature on 
static pressure. It can be seen that 
combustion is possible down to 
temperatures of about 700 K at a pressure 
of 250 kPa. As the pressure drops the 
minimum temperature required to support 
combustion increases. At 50 kPa 
combustion is not possible at 
temperatures below 1200 K. 

It is also noticeable that when 
conditions are such that combustion can 
occur at temperatures below 1000 K, very 
high values of specific impulse are 
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achieved. This is thought to b e due to 
the low coebustion teepersture which 
suppresses dissocistion of the combustion 
products . 

Coabustion could be induced at 
teaperstures below 1000 K either by aesne 
of high pressures or by the use of a Ion# 
coabustion chaaber, and in both cases 
high values of specific iapulse were 
achieved. Fig 6 shows previously 
published data from Ref 7. for an 
intake pressure of 150 kPa at Hach 3. S. 

It is seen that the longer duct produced 
siailar results to the short duct at a 
higher pressure. 

In Fig 5 it can be seen that in the 
region between 700 K and 1000 K the 
aaount of coabustion , as indicated by 
specific iapulse. is very pressure 
sensitive. This is possibly because 
coabustion is aarginal near the ignition 
limits, and a saall rise in pressure can 
markedly reduce the Ignition delay tiaes 
and perait sore heat release within the 
duct. However, if the conditions are such 
that substantial coabustion occurs, that 
is for long ducts or high pressures, then j 
the strong pressure dependency 
disappears. This say be seen by coeparing 
the peak values of specific iapulse in 
Figs 5, 6 and 7. 

In Fig 6 the effect of extra coabustion 
chaaber length is shown for the high 
pressures produced by the 5 degree scoop. 
It is seen that the only difference is a 
slight lowering of the ignition 
teaperature, and very high levels of 
specific iapulse at the teaperatures Just 
above quenching. The data aay be seen to 
be all following the saae pattern, with 
the pressure and coabustion chaaber 
length determining the point of departure 
froa the coaaon curve. The shape of the 
curve is siailar to the idealised 
condition of equilibriua coabustion also 
Shown in Fig 6 

At higher enthalpies. corresponding to 
intake teaperatures above 1500 K, the 
perforaance appears to be fairly 
independent of pressure. An exception to 
this is given by the 50 Kpa intake 
pressure condition shown in Fig 5. which 
shows signs of reduced output at 

teaperatures up to 1800 K. 

In Fig 9 three coabustion zones are 

loosely defined in the pressure 
teaperature plane, using the specific 
iapulse results. In region 1 no 
coabustion is possible for the 

configuration chosen In region 2 

coabustion occurs. but with a strong 


pressure dependence due to ignition delay 
effects. Region 3 includes those cases 
where coabustion is pressure independent 
The exact boundary between regions 2 and 
3 at higher teaperatures la unclear, but 
it appears to lie between 50 and 75 Kpa 
Figure 9 refers only to the short 
coabustion chaaber . There is not enough 
data to fully define a siailar curve for 
the longer coabustion chaaber, but the 
available data suggests s shift of all 
the boundaries to lower teaperatures 
occurs . 

It should be stressed that Fig 9 only 
applies to the specific experiaentai 
configuration used. Its aain value lies 
in providing data to teat theoretical 
aodels. rather than being representative 
of a flight situation. It does indicate 
however, that terse Jets designed to 
optiaiae specific iapulse are likely to 
be operating in pressure sensitive 
regions. It also suggests that the low 
pressures sssocisted with very high 
altitude flight aight experience ignition 
problems, even at high temperatures. 


HUM R I CAL SIMULATION. 

An atteapt was Bade to aodel the mixing 
and coabustion processes in the scraajet 
using s two diaensional parabolic 
computer code. The program is baaed on 
ref 3. but incorporates finite rate 
cheaistry. It is two diaensional in that 
it considers transverse gradients of 
teaperature. cheaical species 

turbulence end velocity. However, it 
cannot aodel • transverse pressure 
gradients, which have been shown in Ref 1 
to be significant in the thrust producing 
processes of two diaensional acraajeta. 

In Fig 10 a acheaatic la shown of the way 
in which the coabustion zone is expanded 
in order to develop thrust. Most of the 
coabustion occurs in a liaited aixing 
region around the centreline, upstrean of 
the expansion propagating froa the start 
of the diverging section. In this narrow 
ducted region it aay seea not 
unraasonabla to aaauae that the effects 
of coabustion are felt uniformly across 
the section of the duct, and that a one 
diaensional pressure treatment aay give 
sufficient eccuracy . 

The authors, in as yet unpublished work, 
have developed a procedure for Batching 
the output froa the quasi two dimensional 
mixing program with linearised small 

perturbation wave theory in the expanding 
region. This has produced good 
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correlations both with the experiment ally 
recorded pressure profiles, and with the 
calculated values of specific iepulse. 

A different approach was adopted for this 
paper, due to the large amount of data to 
be processed. By inspection of Fig 10 it 
can be seen that the exhaust nozzle flow 
will be expanded by multiple reflections 
of the 15 degree corner expansion between 
cowl and the thrust surface. However, 
it was shown in Ref i that valid 
comparisons with a one dimensional code 
can only be made after the elimination of 
all transverse disturbances, that is. in 
the presence of parallel exhaust field. 
For the present configuration this means 
allowing the full reflection of the 
corner expansion from the cowl to 
propagate through the jet. 


The program does predict the sharp drop 
in specific impulse which occurs near the 
ignition point. and the high levels of 
specific impulse observed for the long 
duct at lower temperatures. However, a 
consistently high ignition temperature is 
given, except for the 50 kpa case. The 
high values of specific impulse observed 
•bove the ignition temperature for 
the short duct in the 200 to 500 kpa 
range. Fig 7, were not predicted. 

The results presented in Fig u were 
obtained with oxygen dissociation 
levels as obtained from the non 
•Quilibrim nozzle expansion program. 
Ref 6. Free stream turbulence was 
input as a constant fraction of the 
second power of velocity. 


Therefore the mixing program was run with 
an expanding section such that the area 
ratio experienced by the duct is the same 
as that produced by the double pass of a 
15 degree expansion fan traversing the 
jet. This gives an area ratio of 10.3 for 
the nominal nozzle exit each number of 
3.5. Corrections to this value have not 
been made to account for Mach number i 
variations between conditions. i 

i 

To give a direct comparison with the 
computations . the experimental results 
have to be modified as outlined in Ref i 
to compute the thrust that would have 
realised if all the disturbances 
eliminated by contouring of the 
thrust surface. This has not been done in 
this instance, because the technique of 
Ref l does not apply to long combustion 
chambers with a significant pressure rise 
upstream of the corner. The data is 
♦ therefore presented not for direct 
quantitative comparison with experiment, 
but rather to indicate trends, and to 
predict theoretical maxima for a perfect 
expansion without wave drag. 

In Fig li the results of the computations 

are summarised. 

Quantitative trends which are evident in 
the experimental data also appear in 
the computations. Thus, at a given 
static pressure. the long duct yields 
higher specific impulse than the short 
duct at temperatures just above ignition, 
but at higher temperatures the length of 
the duct does not make much difference 
Also, with a given duct length, the 
*P* c ifl c impulse falls as the pressure 
falls . 


Agreement with individual data points can 
always be obtained by adjustment of 
various parameters in the program, such 
•* free stream turbulence and oxygen 
dissociation levels. The effect of the 
boundary layer which forms on the outside 
walls of the injector strut is also 
unknown, but viscous heating may be 
important in determining the free stream 
ignition temperature. 


only covers that section of the jet which 
is expanded by the first pass of the 
expansion fan, however as explained above 
the computations had to be performed on a 
double pass in order to produce a 
^* r *3^**’ ***** outlet. The computations 
would therefore be expected to 
overpredict the thrust, because of the 
jet expansion which occurred in 
uninstrumented regions of the model 
.which did not contribute to the measured 
thrust production. 


Not withstanding the prediction of 
qualitative trends, it is clear that the 
computational model fails in that it 
underpredicts specific impulse, and does 
not give a correct ignition temperature 
There are two possible reasons for this 


The first is that the 
chemistry models used in the 
not be adequate, and as 
ignition point is concerned 
most likely cause 


mixing and 

prograa may 

fsr as the 
this is the 


The second is that while the one 
dimensional approximation may give 
reasonable results in the constant area 
duct section, and the unpublished results 
mentioned above suggest that this is so 
s one dimensional treatment may not be 
correctly used to model the expansion of 
the expansion traverses the 
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let different point* at the **■• axial 
location will aee different pre.iure., 
and thia eight have a aignificant effect 
on thruat production. 


COWCLUSIQWS, 

Scraajet performance ahow* a atrong 
preeaure dependency in region* near to 
the ignition lieit*. In region* of 
aaxiaua thruat. which are deairable 
target operating condition*, preaaure 
effect, are .1.0 evident at the lower 
te.peraturea At higher enthalpie. ie 
at intake temperatures above 1500 K. 
combustion is insensitive to preaaure 
above intake pressures of 50 Kpa 


Nuserical modeling procedures are not yet 
adequate for quantitative thruat 
predictions , but can give 
indication of general 


a useful 


trenda . 
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Pressure Scaling in the Scramjet Model 


1. Introduction 

Previous experimental data [1] obtained in the shock tunnel 
shows a range of intake temperatures where the integrated pressure 
profile (specific impulse ■ I p ) is very pressure sensitive. For 
higher temperatures, I p was insensitive to pressure, while for 
lower temperatures, the ignition limits were pressure sensitive. 
These results may be interpreted in terms of the reaction and 
ignition time correlations with pressure and temperature reported 
in [2] and [3]. The correlations given in [3] are 


II 

8x10’** exp(9600/T) 

(i) 

V 

P 


Is . 

0.000105 exp( -1 . 12T/1000 ) 
— - 

(2) 

V 

p 



where t i = ignition time (sec) ^ * ignition length (m) 

x - reaction time (sec) 4 R = reaction length (m) 

R 

V = flow velocity (m/s) 

T = inlet temperature (K) 
p = static pressure (atm) 

Hence, near the ignition limits, where is close to the length 

of the model l , a decrease in p will increase and result in no 
m 

combustion within the model. Also, a decrease in T will lead to 

an (exponential) increase in and a very sharp drop in I p . 

In the range of conditions where the fuel always ignites 
within the model length, I p will be most sensitive to pressure 
when 4 R is greater than Hence, above a certain temperature 

(found to be approximately 1400K in [1]), I p will be relatively 
insensitive to changes in p. 


The aim of this experiment (to be performed in August 1988) 
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“ t0 ' XP ' ri " nt * Uy th... pressure-length correl.tlnn. 

U), (2) in th. T4 shock tunnel, with the relatively high st.tic 
pressure, ..berated bp the shock tunnel. .. hop. to very the duct 
pressure in th. combustor .od,l by a lector of ,0 snd be eble to 
observe a pressure rise due to combustion. 

2« The E xperiment- 

The experiment will involve the ...surement of pressure 

distributions in three constant-are. combustion chamber models 

Kith different intake pressures but otherwise eguivalent flow 

conditions (i.e. .... .. ch nu „ ber M , te.per.tur. T and species 

concentrations,. Each model will be scaled geometrically to give 

P^ m constant where *,*1.0 for the largest model (also the lowest 
pressure ) • 

There are. however, problem, involved in getting th. shock 
tunnel to provide ..bivalent flows with a rang. „f pr.s.ur.,, ind 
there are also problems i„ manufacturing three (geometrically 
similar model, with an order- of- magnitude difference in site. 

At. on. end, th. largest model ha. to fit into th. ev.il. hi. test 
facility and be short enough to establish stead, flow within the 
test time. At the small end, we cannot manufacture 
injector struts smaller than 3-4 mm thick and have them operate 
successfully at th. high fuel injection pressures retired. 
Fortunately, (i„ the light of comments made on source flow nossles 

U)> w. expect that th. experiment will be insensitive to ...n 
differences in injector scale. 


The more difficult problem concerns th. generation of 
different pressure flows. Inlet flow for th. small , od .i (hl , h p , 
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will be provided by expanding the test gas from the shock tube 
reservoir through a mach 4 nozzle* It will then directly enter 
the model intake as shown in figure 1. The larger models require 
the same parallel flow but a lower pressures. We considered three 
proposals for generating these low pressure flows... 

1. Run the shock tube at different reservoir conditions and 
expand the flow through the same mach 4 nozzle. Unfortunately, 
changing the pressure in the shock tube reservoir while keeping 
the temperature constant leads to a change in entropy (see e.g. 
the Mollier diagram for air in chemical equilibrium [5]). Harris 
& Warren' [6] have shown that the gas composition is strongly 
correlated with the entropy in the reservoir and, for a factor of 
10 change in pressure, we feel that the corresponding change in 
free radical concentrations is too great for this experiment. 

2. Run the shock tube at the same nominal conditions (as for the 
high pressure case), expand the gas through a first nozzle, 
collect it in a second stagnation region (so that there is a drop 
in total pressure through the normal shock preceding this region) 
and then expand the gas through a second contoured nozzle to the 
final flow conditions. This proposal suffers a similar problem to 
proposal 1 as the gas (with increased entropy) will reach chemical 
equilibrium in this second region. 

3. Run the shock tube at the same nominal conditions, expand the 
gas to a higher mach number (say 8-10) and pass the chemically 
frozen gas through a set of compression wedges to give M=4 at the 
model intake. The shocks from the compression wedges provide a 
drop in total pressure (and increase in entropy) and hence a drop 
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in static pressure relative to the ...11 ..del uithout . chang . 4|> 

the ... composition there sr. relatively low te.p.r.ture. 

behind the oblige shocks). The dis.dv.nts*. her. is thst. to 

provide adequate tot.l pressure drop .cross the compression 

shocks, we initially need to expend the |.s to high msch numbers 

(approx. 10) or, for .oderste. values of H (spprox. 8) we need to 

"shock down” to low ..Ch numbers st the model inlet. The former 

option requires extremely l.r*. nozzle .»d correspondingly 

l.rge set of compression wedges (which sre too large for the 

available facility) to provide adequate flow for the l.rge 
model intake. 

The final configuration for the large model (low pressure) is 
shown i„ figure 2 . y, run the shoc „ ^ ^ ^ ^ 

condition, as for the high pressure model, expand the test gas to 
K*8, compress it to M=3 before the model intake, and then expand 
it to M*4 within the .od.l (but before the point of fuel 

injection). This provides us with a design that will , Ju .t) fit 
within the current facility. 


The operating condition for the small . od el will have 

(approximately , T=1600K, ve3 2 86m/s. , pel.l7.tm and . 

dissociation fraction for oxygen a^i.ee*. These valu „ „ er<> 

computed using a one-dimensional code (7). Table 1 provides a 

more complete list of the test conditions At .k- * 

unions. At this temperature 

•nd dissociation level, w. expect ignition to occur at (or very 

clcse to) the injector. Although the flow . m ba thr „_ 

dimensional and involve turbulent mixing of the fuel jet, w. 

that (if the mixing effects scale with the Reynolds number) 
-c Will be able to decouple the mixing effect, fro. the chemical 
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kinetic effects- by having the same Reynolds number in each model. 

The ignition and reaction lengths (as measured by temperature 
rise) were computed for a premixed flow using a one-dimensional 
chemical kinetic code [8]. For a fuel equivalence ratio *=1.0 and 
the flow conditions mentioned above , the numerical code gives 
l (p=l . 1 7 ) =0 . 095m and l (p=0 . 2 ) =1 . 9m. These values may be 
compared with 4 R < p=l . 1 7 ) =0 . 044m and * R < p=0 . 2 ) =0 . 887m computed from 
the correlation (2). Although the lengths are different by a 
factor of 2, the equivalent exponent of pressure is only IX 
different. The numerical simulations also confirmed that the 
reaction length was insensitive to (see also ll]). 

Measurements of the reaction length in the models will be 
obtained indirectly by measuring the static pressure at the wall 
at several downstream locations. To obtain sufficient spatial 
resolution we may have to construct a single pressure profile from 
several shots with the pressure transducers at slightly different 
locations for each shot. However, we will not be able to avoid 
the confusing pressure disturbances introduced at the trailing 
edge of the injector and caused by the mismatch in pressure of the 
fuel jet and the free stream. If it can be arranged, we might 
also try to measure ignition lengths from optical emissions. 

3. Results 

. . . should be interesting when we get them. 


4. References 


1. R.G. Morgan &. R.J. Stalker 1987 : "Pressure scaling effects in 
a scramjet combustion chamber". 8th Int. Symp. on Air Breathing 
Engines, Cincinnati Ohio, June 15-19. 
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Table 1 : Hypothetical Test Conditions 


Paraaeter 

Units 

Sb&II 

noule before 
exit injection 

Model 
Med i ua 

nozzle before 
exit injection 

Large 

nozzle before 
exit injection 

M 


4.23 

4.17 

5.64 

4.17 

8.16 

4.17 

p, static 

ftts 

1.17 

1.17 

0.174 

0.4722 

0.0138 

0.114 

T 

K 

1599 

1630 

1001 

1620(0.6) 

501 

1683(3) 

p xlO* 

g/cc 

2.56 

2.51 

0.612 

1.026 

0.0965 

0.251 

V 

l/s 1 

3286 

3274 

3500 

3295(0.6) 

3649 

3326(2) 

a 

X 

1.68 

2.02 

2.67 

o 

7 

— 

1.31 

1.34 

1.39 

l 

■ 

0.1 

0.4 

l.C 

1 

■ 

6 

w 

deg 

— 

• 

12. 

.0 

15. 

0 

e. 

i 

deg 

10.0 

10. 

.0 

10. 

0 


Quantities in parentheses are percentage differences between the current aodel 
and the snail aodel. 

Shock tube stagnation conditions are H q = 7.5 MJ/kg, P Q = 450 ata, and 

T = 5000 K. 
o 
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Abstract 

This pa par praaa&ts the preliminary test flow 
raaulta froa the University of Queensland's fraa 
piston expansion tuba. It has baan shown that pitot 
pressure naasuranants ara qualitatively sinilar to 
thosa Iron tba Langley facility. Furthernore, 
additional tast conditions which bars acceptably 
steady flows of air are demonstrated to exist using 
pitot prassura naasuranants. Shortened test times 
are explained partially by contamination from low 
density driwer gas bubbles and reflection of 
expansion wares iron tba driver-test gas interface. 


Introduction 

A renewed interest in hypersonic flight bas 
naturally created a renewed interest in test 
facilities which ara capable of producing high 
enthalpy flows. Tba expansion tuba is an impulse 
facility which in theory is capable of sinulating a 
wide range of hypersonic conditions. Bowever, an 
extensive investigation made by Hiller (1917) on an 
At^Fiven expansion tuba found only a narrow band 
of conditions for each test gas in which the tast 
flows were acceptably steady. At the University of 
Queensland an expansion tuba was built with a free 
piston driver in order to determine it the 
versatility of the free piston driver would allow 
the 'window of test conditions' to he widened. 

This paper reports on experiaents performed 
with the University of Queensland expansion tuba. In 
these experiaents the objectives were 

(i) to obtain a coaparsion of the centre-line 
pitot pressure aeasureaents of the 
University of Queensland and the Langley 
tuba for similar running conditions , 

(ii) to show that different test conditions 
could be obtained using different driver 
gases , 

(ill) to show test conditions of the same 
stagnation enthalpy could be obtained 
with different test gas pressures. 


Description of Facilities 

Langley's Six Inch Ixpansion Tube 

A full description of the Langley six inch 
expansion tube is given by Hoore (1975). Briefly 
the driver is jt^Ha in length with an internal 
diameter of 0.356 n. The shock tube or driven 
section is 4»ti a long with an inside diaaeter of 
0.152 n. The acceleration tube is iyn long with 
the same inside diameter. Both the shock tube and 
acceleration tube could be varied in length. 3 The 
duap tank has an approximate volume of 12.8 a and 
is 10.7 n long with an internal diaaeter of 1.22 m. 

♦ Research Associate. 

* Prof. Space Kngineering. 

I Research Fellow. 


Umi varsity of Queensland's Sxpansioa Tube 

The Qniversity of <*eenoland expansion tube 
(see figure 1) has a free piston driver which is 

2.3 a long, 0.1 m internal diameter and uses a 

3.4 kg piston. The shock tube and acceleration tube 
are 2.08 n and 3.18 m long respectively amd both 
have an internal diameter of 0.0386 a. The 
acceleration tube empties into a dump tank-te^t 
section with a volume of approximately 0.15 n. 
Initially the driver and shock tube are separated by 
the primary diaphran and the shock tube and 
acceleration tube are separated by the less 
substantial secondary diaphran. The primary 
diaphran was either 0.1 am thick cold roled mild 
steel or a. 57 mm thick aluminium. The secondary 
diaphram was 0.023 an thick cellophane. The 
acceleration tube and dump tank are reduced to the 
sane initial pressure. 

The shock speeds in both the shock tube and 
acceleration tube were measured using ionization 
gauges. Three gauges were used in the shock tube 
and four were used in the acceleration tube. Froa 
these gauges both the speed and attenuation of the 
shock can be measured. 

Two pressure traces were recorded every run. 
The wall pressure in the shock tube immediately 
upstrean of the secondary diaphran was observed in 
addition to the pitot pressure at the exit of the 
acceleration tube. Bach record consisted of 2048 
data points and was sampled at a rate of one 
megahertz. 


general Theory 

Figure 2 displays the wave diagran for 
operation of an expansion tube. Viscous effects 
have been neglected. It shows that a shock produced 
by bursting the prinary diaphran propagates at a 
constant velocity through the shock tube until 
striking the secondary diaphram. This diaphran then 
bursts and the shock is instantaneously accelerated 
and subsequently moves at a constant velocity down 
the acceleration tube. It is followed by the 
test -acceleration gas interface and the unsteady 
expansion produced from the bursting of the 
secondary diaphram. Hiller (1977) notes that the 
naxinun theoretical teat time is that time between 
the arrival of the test gas and the unsteady 
expansion. BoweVer, this can be decreased. 

The test time can be shortened by the existence 
of "bubbles" or regions which are lower in density 
than the test gas (gourlay (1988)). Their presence 
results in large oscillations of the pitot pressure 
and their arrival makes the test flow useless. It 
is believed that these bubbles result froa mixing at 
the driver-test gas interface and are regions in 
which the driver amd tost gases continually mix. 
When these low density bubbles encounter the 
unsteady expansion they accelerate through its 
pressure gradient and emerge into the test gas with 
a higher velocity than both the velocities of the 
unsteady expansion and the test gas. Bence, it the 
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(i) Figure * displays the pitot pressure when 
the shock tub# was filled to 3.5 kPa ud the 
acceleration tub# was filled to 11 Fa. t helium 
driver was us#d. Fro# tb# shock sp##d is tb# 
acctlaratiom t#b# It can b# shown that total 
enthalpy of this tast flow is 42 MJ/kg. It can b# 
s#aa that a wary steady flow with a total pressure 
of 100 kFa exists tor 10 m§, aft#r which bubbles of 
driver gas and tb# aff#cts of tb# expansion appear. 

(11) bn air driw#r gas was us#d to produc# tb# 
pitot pressure tract giw#n in figur# 4. Tb# shock 
tub# filling pressure was 1.1 kPa and tb# 

acceleration tube was filled to 1.3 Pa. The total 
enthalpy was 11 MJ/kg. It can be se#n that a 
reasonably steady flow exists for approxinat«ly 
200 ms for a total pressure of 128 kPa. 

(iii) Figure 1 displays the b#st pitot pressure 
n#asur#eent produced with an argon driwer . The 
shock tube filling pressure was 3.5 kPa and the 
acceleration tube was filled to 33 Pa. Although 
this result is not as clean as that for the helium 
driwer it can be seen that after the initial rise 
due to the acceleration gas there is approximately 
100 Mi of test tine in which the pitot pressure 
vould wary no nore than ±5 % f roe 535 kPa. The 
total enthalpy for this test condition is 8.4 MJ/kg. 

These conditions were repeatable approxinately 
three out of ewery four tines. 

ZZZ lew Test Flows bt Different Pressures 

If the test tine is indeed limited by the 
mechanism outlined above, then the characteristics 
in figure 2 would be the sane provided the ratios of 
(i) the primary diaphran burst pressure to the 
shock tube filling pressure, and 
(11) the shock tube filling pressure to the 
acceleration tube tilling pressure 
are presarwed. Thus, in theory, similar gualitatiwe 
pitot pressure traces should be seen if these ratios 
remain static. Only the pressure level should 
change. The shock speeds' and thus the total 
enthalpy should remain the sane. 

This conclusion was tested by performing 
experiments at reduced operating pressure levels. 
Figures 8e-8f are the pitot pressure traces when the 
acceleration tube filling pressure was varied from 
2.7 Fa to 133 Fa and the shock tube filling pressure 
was 1.7 kFa. The primary diaphran bursting pressure 
was approxinately 11.5 MPa and the driver gas was 
argon. Zn figures 9a-9f the shock tube filliog 
pressure was 3.4 kFa and the acceleration tube 
filling pressure was varied between 5.3 kPa and 


246 kFa. The primary diaphran bursting pressure was 
35 KFa and the driwer gas was again argon. The 
above ratios ((1) end (ii)> have been basically 
preserved but the overall pressures have doubled. 

Zt can be seen from these preliminary results 
that the inviscid theory would appear to be correct, 
bt the lowest acceleration tube filling pressure 
there is a gradual rise in pressure followed by a 
relatively noisy flow that is believed to result 
from the lateral pressure oscillations produced by 
bursting of the secondary diaphran. hs the 
acceleration tube filling pressure is Increased the 
reflection of the expansion arrives earlier and the 
step due to the acceleration gas becomes more 
pronounced. Furthermore the acceptable test time 
increases from figures Sa and 8a to 8c and 9c and 
then decreases as the effect of the bubbles and 
reflected expansion become more and more dominant. 
The pressures of the plateaus of figures 8c and 9c 
are 140 kFa and 280 kFa respectively. The 
enthalpies of these flows are 8.4 MJ/kg and 
8.2 MJ/kg respectively. 


Conclusions 

(i) The Langley facility and the University of 
Queensland facility are qualitatively comparable for 
similar teynolds numbers and shock tube shock 

speeds. 

(ii) Different driwer gases can be used to 
obtain steady flows for different stagnation 
enthalpies. 

(iii) Preliminary results indicate that 
different operating pressure levels can be used for 
similar stagnation enthalpies provided the ratios of 
the diaphran burst pressure to the shock tube 
filling pressure and the shock tube tilling pressure 
to the acceleration tube filling pressure are 
maintained constant. 
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riGURE 7. Pitot Pressure Trace for an Argon Driver with P^ - 3.5 *Pa and P fc • 33 Pa 
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FIGURE 8. Pitot Pressure Traces. Ar Driver 

P > 1.7 kPa 
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FIGURE 9 


Pitot Pressure Traces. Ar Driver 

P - 3.4 kPa 
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